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5 CROSS REFERENCE TO APPENDICES 

This application is related to Application Attorney Docket No. 10001606-1 
which is incorporated herein by reference in its entirety. 

CROSS REFERENCE TO APPENDICES 

Appendix Al 5 which is part of the present disclosure, is an appendix 

1 0 consisting of 8 pages. Appendix Al lists source code written in C programming 
language of an illustrative embodiment of the present invention using an eight-bit 
code to control a printer system having a pulse width modulated laser capability for 
each printed pixel. Appendix A2, which is part of the present disclosure, is an 
appendix consisting of 13 pages. Appendix A2 lists the particular coefficients used in 

1 5 the Tone-Dependent Error Diffusion process for the eight-bit code used in accordance 
with an embodiment of the present invention. Appendix Bl 5 which is part of the 
present disclosure, is an appendix consisting of 14 pages. Appendix Bl lists source 
code written in C programming language of an illustrative embodiment of the present 
invention using a two-bit code to control a printer system having a pulse width 

20 modulated laser capability for each printed pixel. Appendix B2, which is part of the 
present disclosure, is an appendix consisting of 2 pages. Appendix B2 lists source 
code written in C programming language that is used in conjunction with Appendix 
Bl for a two-bit code used to control a printer system having a pulse width modulated 
laser capability for each printed pixel. Appendix B3, which is part of the present 

25 disclosure, is an appendix consisting of 22 pages. Appendix B3 lists the particular 
coefficients used in the Tone-Dependent Error Diffusion process for the two-bit code 
used in accordance with an embodiment of the present invention. 

COPYRIGHT NOTICE 

A portion of the disclosure of this patent document contains material that is 
30 subject to copyright protection. The copyright owner has no objection to the facsimile 
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reproduction by anyone of the patent document or the patent disclosure, as it appears 
in the patent and trademark office patent file or records, but otherwise reserves all 
copyright rights whatsoever. 

FIELD OF THE INVENTION 
5 The present invention relates to halftoning, and in particular to a halftoning 

method that combines dot size modulation and dot spacing modulation to control the 
overall gray level in an image forming or printing device. 

BACKGROUND 

Continuous-tone images, such as charts, drawings, and pictures, may be 

10 represented as a two-dimensional matrix of picture elements (pixels). The spatial 

resolution and intensity level for each pixel are chosen to correspond to the particular 
output device used. Typically, digital halftoning is used to transform a continuous- 
tone image into the desired matrix of pixels. 

Conventional methods for digital halftoning generally fall into two categories: 

15 clustered dot and dispersed dot. As is well known in the art, in clustered dot 

techniques, the size of the printed dot is varied to control the perceived gray level or, 
equivalently, density of the printed tone. These methods may be thought of as 
amplitude modulation (AM) halftoning techniques since the amplitude or size of the 
dots controls the printed gray level. 

20 Fig. 1 is a block diagram illustrating a conventional clustered dot halftoning 

process. As shown in Fig. 1, typically, an image x(m, n) is transformed, or tone 
compensated, at block 2, so that the printed gray level for each input value is correct. 
Tone compensation (block 2) is used to correct for attributes of the halftoning 
algorithm, image calibration, and printer behavior. In particular, most printers are not 

25 linear due to the inevitable effects of dot overlap. The dot size is then modulated at 
block 4 and the resulting signal is sent to the printing system. Typically, dot size 
modulation is performed using screening, which uses a threshold array that is 
compared to the pixel value. AM halftoning methods have certain advantages such as 
increased robustness to printing artifacts. However, a limitation of AM halftoning 

30 methods is that the choice of dot size presents a fundamental tradeoff between spatial 
resolution and number of gray levels. 
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Dispersed dot halftoning methods control the printed gray level through the 
spacing or, equivalently, the frequency of dot placement. Dispersed dot halftoning 
may be thought of as a frequency modulation (FM) halftoning technique since the 
frequency or spacing of the dots controls the printed gray level. Dispersed dot 

5 halftoning methods are well known, and include, for example, error diffusion 

halftoning, screening and, most recently, iterative search based halftoning. Fig. 2 is a 
block diagram illustrating a conventional dispersed dot halftoning process. As with 
the clustered dot halftoning process, the image x(m, n) is transformed, or tone 
compensated at block 6, so that the printed gray level for each input value is correct. 

10 The spacing of the dots is then appropriately dispersed at block 8 and the resulting 
signal is sent to the printing system. 

Typically, dispersed dot halftoning uses the smallest dots possible to print as it 
is undesirable to actually notice the dots in the printed image. However, a drawback 
of dispersed dot halftoning is that the printer must be capable of printing well formed 

15 isolated dots. For example, error diffusion is typically used only in ink jet type 
printers because they generally can print stable isolated dots. However, error 
diffusion is typically not used in commercial electrophotographic printers and copiers, 
such as laser printers, due to their instability in producing binary or multilevel 
halftones. 

20 SUMMARY 

A halftoning system, in accordance with an embodiment of the present 
invention, uses dispersed dot halftoning in conjunction with dot size modulation to 
produce an image in which the density and size of the dots are modulated in 
conjunction to control the overall gray level in an image forming system. The dot 

25 density is modulated to control the spacing or frequency of the dot to be printed for a 
particular pixel location with respect to dots to be printed at preceding and subsequent 
pixel locations, while the dot size is modulated to control the size of the dot to be 
printed at the particular pixel location. The dot density and size modulation system 
may be used with an image forming or printing system, e.g., a computer and a printer. 

30 It should be understood that the image forming or printing system of the present 

invention is intended to include, but is not limited to, an electrophotographic printer 
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or copier, an ink jet printer, a facsimile machine or any other device that may be used 
to print an image. The dot density and size modulation system is particularly useful in 
modern electrophotographic printing systems that allow the printed dot size to be 
almost continuously varied through the specification of a pulse width modulation 
5 (PWM) code. 

In accordance with an embodiment of the present invention, the dot density is 
modulated by controlling the dot density for a pixel location using an input pixel 
value for that pixel location and by performing dispersed dot halftoning to produce a 
dot position based on the dot density. For example, a look up table may generate a 

10 dot density value based on the input pixel value, which is then used in the dispersed 
dot halftoning process. The dispersed dot halftoning process may be error diffusion, 
e.g., tone dependent error diffusion, dispersed dot screening, or iterative search based 
halftoning. Tone dependent error diffusion is particularly advantageous because it 
generates high quality dispersed dot patterns and is computationally efficient. 

15 The dot size is modulated, in accordance with an embodiment of the present 

invention, by independently controlling the dot size for the pixel location using the 
input pixel value and performing dot size modulation based on both the dot size and 
the dot position. Thus, for example, a separate look up table may be used to generate 
the dot size value based on the input pixel value. Advantageously, the dot size control 

20 and dot density control may be precomputed, and in one embodiment are designed 
with respect to each other and the characteristics of the printing system to optimize 
the image while matching the printing device characteristics. The resulting dot size 
value, as well as the dot position from the dispersed dot halftoning process, are then 
used to generate a signal, e.g., a pulse width modulated code, that controls the printer 

25 system. The signal may be produced by another look up table and may include 

information regarding both the pulse width of the desired dot for the pixel location as 
well as the justification, e.g. left, center, or right, for the pixel location. 

In accordance with an embodiment of the present invention, the dot size look- 
up table and dot density look-up table are designed with respect to each other and the 

30 characteristics of the printing system in a closed loop measurement process to 

optimize the quality of the image while matching the printing device characteristics. 
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The system, in accordance with the present invention, uses a combination of 
the dot size and dot density to produce the desired gray level, which offers advantages 
over pure dot density modulation systems or pure dot size modulation systems. The 
combination of dot density and dot size modulation allows an extra degree of 
5 flexibility, which can be used to increase the visual quality of the halftoned pattern 
and/or increase the robustness of the halftoning to printer artifacts and variations. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustrating a conventional clustered dot halftoning 
process. 

10 Fig. 2 is a block diagram illustrating a conventional dispersed dot halftoning 

process. 

Fig. 3 is a block diagram illustrating a dot density and dot size modulation 
(referred to herein as AM/FM halftoning) system, in accordance with an embodiment 
of the present invention. 
1 5 Fig. 4 is a schematic diagram of a printing system, including a printer and a 

computer that may be used with the present invention. 

Fig. 5 shows a block diagram of an embodiment of the AM/FM halftoning 
system in which the dot density control and dot size control are look-up tables and the 
dispersed dot halftoning is tone dependent error diffusion. 
20 Fig. 6 shows a block diagram of a tone dependent error diffusion, which may 

be used in the AM/FM halftoning system of Fig. 5. 

Figs. 7A and 7B show an array of pixels and a pixel pair, respectively, 
illustrating the error diffusion used in accordance with an embodiment of the present 
invention. 

25 Figs. 8A, 8B, and 8C show an idealized laser modulation output for a PWM 

system in which the width of the pulse may be less than the full width of the pixel and 

may be left, center, or right justified, respectively. 

Fig. 9 illustrates the method used to form a circular dot with an area 9=2 and 

that is centered at pixel (m,n). 
30 Fig. 10 is a flow chart showing a method of optimizing the dot size and dot 

density look-up tables using a closed loop measurement process in accordance with an 

embodiment of the present invention. 
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Fig. 1 1 is a graph showing the measured tone curve, ge(n), i.e., the output 
absorptance, versus input dot density for a fixed value of dot size 0 = 1, where a 
square dot shape is used. 

Fig. 12 is a graph showing the print distortion, De(n) versus input dot density 
5 for a fixed value of dot size 0=1, where a square dot shape is used. 

Fig. 13 is a graph showing a set of inverse tone curves for dot sizes ranging 
from 0.3 to 1.4. 

Fig. 14 is a graph showing a set of plots of the print distortion versus the 
desired output absorptance for dot sizes ranging from 0.3 to 1.4, where each curve 
10 corresponds to De jg as a function of g for fixed 6. 

Fig. 1 5 is a graph showing a dot size curve for different a values, where the 
optimized dot size look-up table, 0i used is a = 0.05. 

Fig. 16 is a graph showing a dot density curve for an optimized dot density 
look-up table that corresponds to the curve a = 0.05 shown in graph in Fig. 15 and is a 
1 5 function of the desired absorptance. 

Fig. 17 is a graph showing another dot size curve for an optimized dot size 
look-up table as a function of desired absorptance. 

Fig. 18 is a graph showing a dot density curve for an optimized dot density 
look-up table that corresponds to the curve shown in graph in Fig. 17 and is a function 
20 of the desired absorptance. 

Fig. 19 shows another embodiment of the AM/FM halftoning system, in 
accordance with the present invention. 

Fig. 20 shows a block diagram illustrating another embodiment of the AM/FM 
halftoning system. 

25 DETAILED DESCRIPTION 

Fig. 3 is a block diagram illustrating a dot size and density modulation system 
(referred to herein as an AM/FM halftoning system) 100, in accordance with an 
embodiment of the present invention. As shown in Fig. 3, each pixel x(m,n) of an 
image to be printed is processed to determine the dot density, at dot density control 

30 102, and the dot size, at dot size control 104. The AM/FM halftoning system 100 
includes a subsystem 106 with a FM modulation portion, , e.g., dispersed dot 
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halftoning unit 108, and an AM modulation portion, e.g., dot size modulation unit 
110. 

The dot density control 102 provides a dot density value that is received by the 
dispersed dot halftoning unit 108 to determine the positions of dots on the discrete 
printing grid. The dot size control 104 provides a dot size value that is received by 
the dot size modulation unit 110 along with the positions of the dots from the 
dispersed dot halftoning unit 108. The dot size modulation unit 110 generates an 
output signal indicating the modulated size and modulated density of the dot to be 
printed for the input pixel x(m,n) by the printing system. 

Fig. 4 is a schematic diagram of a printing system 112 that includes a host 
computer 1 13, a monitor 1 14 (e.g., a CRT), and a printer 116, which may be used 
with the present invention. Printer 116 could be any type of printer that prints black 
and/or color dots, including an electrophotographic printer, e.g., a laser printer, or 
copier, or an inkjet printer. Printer 116 typically includes a toner cartridge 117, which 
delivers particles of dry toner via electrostatic attraction to a rotating drum or sheet of 
paper. The paper is then heated to melt the toner so that the toner permanently 
adheres to the paper. For inkjet printers, a print cartridge is typically used that scans 
across the paper and prints droplets of ink in a well known fashion. Additional toner 
cartridges 118 may also be included in printer 116, for example, to provide different 
color inks. Printer 116 often includes a printer controller 1 19 for controlling the 
printing of dots, and in fact, may include the host computer 113 within printer 1 16, in 
which case monitor 1 14 is not needed. There are many well-known types of these 
devices, and details of their operation need not be presented here for a clear 
understanding of the operation of the present invention. 

Fig. 5 shows a block diagram of one embodiment of an AM/FM halftoning 
system 120 in accordance with the present invention, in which the dot density control 
is a first 256x7 bit look-up table ("LUT") 122 and dot size control is a second 256x7 
bit LUT 124. In this embodiment, the dot density control LUT 122 and dot size 
control LUT 124 map each input pixel x(m, n) to a specific desired value. In one 
embodiment, LUT 122 produces a dot density value ranging from 0 to 128 and LUT 
124 produces a dot size value ranging from 0 to 63, but these values may vary if 
desired. For example, in some printing systems the dot density value may vary 
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between 0 and 255. Moreover, the dot density control and the dot size control, in 
other embodiments, may be more general and may depend on the local neighborhood 
of the input pixel x(m, n), which might be useful, e.g., for purposes of enhancing the 
quality of text or image edges, or more accurately rendering textures or smooth 
5 gradients. 

Advantageously, both the dot density and dot size are varied independently 
based on attributes of the input image x(m, n). Nevertheless, neither the dot density 
control nor the dot size control independently control the printed tone. Instead, the 
printed tone is controlled through a combination of the two. The dot density and size 

10 are designed with respect to each other and the printing system to optimize the printed 
image while matching the device characteristics of the printing system. 
Consequently, tone correction may be incorporated directly into the combined 
density/size control. 

As shown in Fig. 5, the FM modulation the dispersed dot halftoning unit is a 

15 tone dependent error diffusion (TDED) system 128. The TDED system 128 may be, 
for example, similar to that described in U.S. Patent Serial No.09/307,003, entitled 
"Tone Dependent Error Diffusion" to Pingshan Li and Jan P. Allebach, filed May 7, 
1999, which is incorporated herein by reference. Tone dependent error diffusion is 
advantageous because it generates high quality dispersed dot patterns and is 

20 computationally efficient. The dispersed dot halftoning unit, however, may depend 
on the particular application, i.e., printer system, and may be any one of a wide 
variety of known halftoning methods including error diffusion, dispersed dot 
screening, or iterative search based halftoning such as direct binary search (DBS), 
which is described, for example, in M. Analoui and J. Allebach, "Model-Based 

25 Halftoning Using Direct Binary Search,", Proc. SPIE/IS&T's symp., electr, imag. sci. 
and tech., vol. 1666, pp. 96-108, San Jose, CA, Feb., 1992, and which is incorporated 
herein by reference. 

Typically, the quantity varied in dispersed dot halftoning is the dot spacing, 
which is specified using the number of dots per unit area or dot density. It should be 

30 understood, however, that for some dispersed dot halftoning methods the output of the 
dot density control may be only approximately equal to the true dot density. For more 
information regarding dispersed dot halftoning in general, see for example, J. 
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Mulligan and A. Ahumada, Jr., "Principled Halftoning Based on Models of Human 
Vision," Proc. SPIE/IS&T's symp., electr, imag. sci. and tech., vol. 1666, pp. 109- 
121, San Jose, CA, Feb., 1992; and T. Pappas, and D. Neuhoff, "Least-Squares 
Model-Based Halftoning," Proc. SPIE/IS&T's symp., electr, imag. sci. and tech., vol. 
5 1666, pp. 165-176, San Jose, CA, Feb., 1992, which are incorporated herein by 
reference. 

Fig. 6 shows a block diagram of a TDED system 140, which may be used as 
the TDED system 128 of Fig. 5 to produce a TDED output in the form of a halftone 
value, in accordance with an embodiment of the present invention. The input pixel 

10 value, which has continuous-tone value, e.g., from 0 to 128, is represented as X(m,n), 
where m and n are the pixel locations. The TDED system 140 uses an error diffusion 
filter 142, which acts as a weighting matrix, and a threshold matrix 144. 

The input pixel value is placed in a weight LUT 146, which may be, e.g., a 
3x129x8 bit LUT, and which produces three different weights for the error diffusion 

15 filter 142. The input pixel value is also placed in a threshold weight LUT 148, which 
may be, e.g., a 2x129x8 bit LUT, and which produces two different thresholds for the 
threshold matrix 144. In addition, the pixel location (m,n) is received by a direct 
binary search screen (DBS) 150, the output of which is multiplied by the resulting 
output of threshold weight LUT 148, which is then summed with the other resulting 

20 output of threshold weight LUT 148. The summed result is then provided to the 
threshold matrix 144. 

The error diffusion filter 142 provides an error value from previously 
processed pixels that is summed with the input pixel value X(m,n), to produce a 
modified pixel value u(m,n). The modified pixel value u(m,n) is received by the 

25 threshold matrix 144 and compared to at least one threshold level to produce the 
output halftone value g(m,n). A quantizer error value d(m,n) is produced as the 
difference between the output halftone value g(m,n) and the modified pixel value 
u(m,n). The quantizer error d(m,n) is received by the error diffusion filter 142, and is 
diffused to neighboring, subsequently processed pixel locations. See U.S. Patent 

30 Serial No. 09/307,003, entitled "Tone Dependent Error Diffusion" to Pingshan Li and 
Jan P. Allebach, filed May 7, 1999, which is incorporated herein by reference for a 
discussion of the implementation of the TDED system 140. As one of ordinary skill 
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in the art understands, a TDED system may be optimized for particular printing 
systems and printing characteristics by altering the particular process coefficient. For 
the particular coefficients used in accordance with an embodiment of the present 
invention, see Appendix A2. 
5 Figs. 7 A and 7B show an array of pixels 160 and a pixel pair 164, 

respectively, illustrating the error diffusion of the TDED system 128 in accordance 
with one embodiment of the present invention. As illustrated in Fig. 7 A, the dots are 
printed using a serpentine raster scan as indicated by arrows 161 . The pixels are 
organized in pixel pairs, in which a modulated dot will be fired in only one pixel 

1 0 illustrated with an "X". The pixels which are illustrated with an "O" will not have a 
dot fired. The dots are constrained to fire at positions on a diagonal grid pattern. The 
tone dependent error diffusion uses four weights that are diffused to four adjacent 
pixels, as indicated by arrows 162 and 163. As shown in Fig. 7 A, the error is diffused 
to the next pixel to be processed, regardless of whether a dot will be fired in that pixel 

15 and to the three adjacent pixels in the next row, again regardless of whether a dot will 
be fired in that pixel. As shown in Fig. 7B, which shows a pixel pair 164, in one 
embodiment of the present invention, each pixel in a pixel pair 164 is justified 
towards the other pixel. Thus, the left hand pixel "X" is right justified and the right 
hand pixel "O" is left hand justified. 

20 As shown in Fig. 5, the dot size modulation is a pulse width modulation 

(PWM) system 130 that directly modulates the size of the physical dot, as is well 
known in the art. The dot size modulation step requires very little computation since 
the appropriate PWM codes for each dot size may be precomputed and stored in a 
LUT and is accessed using the halftone value from the TDED system 128 and the dot 

25 size value from the dot size LUT 124. The dot size modulation, however, may also 
depend on the particular application, i.e., printing system, and may be any one of a 
number of known methods which vary the size of a printed dot either by grouping 
clusters of dots together, including clustered dot digital halftoning. 

The PWM system 130 produces a signal, e.g., an eight-bit PWM code, that 

30 controls the size of the printed dot, for example, produced by a laser in an 
electrophotographic printer. As is well understood in the art, in an 
electrophotographic printer, the toner is only deposited when the laser is on. Thus, 
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each six-bit PWM code controls the width of a single pixel. The output signal from 
the PWM system 130 may also include two bits of information regarding the 
justification of the dot, i.e., whether the dot is left, center or right justified within the 
pixel. 

5 The source code in Appendix Al written in C programming language is one 

example of an implementation of the present invention using an eight-bit PWM code 
to control a printer system having a pulse width modulated laser capability for each 
printed pixel. 

In Fig. 7 A each grouping of two adjacent pixels is referred to as a pixel pair 
10 and is shown with an X on the left pixel and a 0 on the right pixel. A single example 
pixel pair is indicated at the bottom of Fig. 7. Each pixel pair is either enabled or 
disable depending on the TDED output value at the X position corresponding to the 
left hand pixel in the pixel pair. If TDED output value indicates that a dot should fire 
at the X position, then the pair is enabled. If the TDED output value does not indicate 
15 that a dot should fire at the X position, then the pair is not enabled. If the pair is not 
enabled, then the PWM codes for both the X position and the 0 position of the pixel 
pair are set to 0 and nothing is printed at either of the two pixels. If the pixel pair is 
enabled, then each of the two pixels are AM modulated separately using dot size 
value from the dot size LUT 124. More specifically, the PWM code for the left pixel, 
20 denoted by XPWM, and the PWM code for the right pixel, denoted by OPWM, are 
computed as follows: 

XPWM -DotSizeLUT(left pixel)»l+RightHand Justification; and equ. 1 

OPWM =(DotSizeLUT(right pixel)+ 1 )»1 +LeftHandJustification. equ. 2 

If a dot is not fired, then XPWM-OPWM-off. Here the values "left pixel" and "right 
25 pixel" denote the values of the continuous tone image to be printed at the 

corresponding left and right locations of the pixel pair. 

Figs. 8A, 8B, and 8C shows an idealized laser modulation output for a PWM 

system in which the width of the pulse may be less than the full width of the pixel and 

may be left, center, or right justified, respectively. Thus, for example, the lower six 
30 bits of a eight-bit PWM code may specify the width of the laser pulse, and the highest 

two bits specify the justification (left, center, or right). The pulse width value ranges 
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from 0 to 63 , where 0 specifies a width 0 pulse, and 63 specifies a pulse that is the fall 
width of the pixel. The PWM codes that may be used to specify the justification for 
each modulated pixel are given in Table 1. 



l sl Bit 


2 nQ Bit 


Value 


Justification 


1 


1 


OxcO 


Right 


0 


1 


0x40 


Left 


0 


0 


0x00 


Center 


1 


0 


0x80 


Undefined 




fable 1 



Fig. 9 illustrates another method to form a dot, in accordance with an 
embodiment of the present invention. Fig. 9 shows a method used to form a circular 
dot 170 with an area 9=2 and that is centered at pixel (m 3 n). First, the radius of the 

desired ideal dot is determined as r = 4q/k . Then, for each row of pixels, the area of 
10 intersection between the row and the ideal dot is calculated. To do this, let k index 
the rows of the image, and let k = 0 for the row containing the pixel (m, n). Then the 
region of intersection is defined by 

S k = {(x,y):x 2 +y 2 <r and| x -k |< 1/2} equ. 3 

where x and y are the distance from the center of the pixel (m, n). The area of 
1 5 intersection is then pk = Area{Sk} . If Pk > 1, then nk = 2L(pk- l)/2j + 1 pixels in the 
l^th row fijUy turned on. For the remaining fractional area yk ^Pk - nk, two pixels 
on either side of the fully turned on pixels are turned on, each with area yk/2. The 
group of pixels are justified so that the dot is more compact, i.e., the left side pixel is 
right justified and the right side pixel is left justified so that the dot is formed by one 
20 pulse. If Pk < 1 , we just turn on one pixel of area Pk with center justification. This is 
illustrated in Fig. 9 as array 172. The PWM codes, with values ranging from 0 to 63 
and an indication of left (L), right (R), and center (C) justification, for each pixel is 
shown in array 174 in Fig. 9. Importantly, the mapping of the dot size 0 to the PWM 
codes can be precomputed and implemented using a LUT in PWM 130 (shown in Fig, 
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5). Precomputing the PWM codes and using a LUT increases the speed of operation. 
In some embodiments, a rectangular dot rather than a circular dot may be used. The 
use of a rectangular dot has the advantage that when the dot area 0 is less than 1 only 
a single pixel is modulated. The use of a rectangular dot makes the dot printing more 
5 stable for dot size 9 > 1 . 

When the dot area is greater than 1, i.e. 9 is greater than 1, dot overlap may 
occur. In this case, the pulse widths are added, and only the first pulse justification is 
retained. More specifically, suppose there are N dots indexed by i = 1,—, N that 
correspond to the same pixel location. Let bi be the justification part (the two most 
10 significant bits) of the PWM code, and let a* be the pulse width part (the remaining six 
bits) of an eight-bit PWM code. Then the PWM code for a pixel is determined by 

N 

C m>R = \ + min(63, £ a t ) equ. 4 

The value of the pulse width is clipped if it exceeds the maximum value of 63. The 
effects of this clipping on the tone curve are removed when the dot density control 
15 LUT 122 (Fig. 5) and dot size control LUT (124) are designed to prevent distortion of 
the tone curve. 

In some cases, 63 pulse width levels of PWM may be available or less than 63 
levels of pulse width modulation may be available, which may cause contouring 
artifacts in the printed images. In either case, contouring happens because a change of 

20 one level in the PWM code is substantial enough to cause a visible difference in gray 
level. To eliminate this problem, the pulse width codes may be dithered. The 
dithering may be done using random numbers, a white noise threshold mask, a 
dispersed dot threshold mask (such as the DBS screen 150 shown in Fig. 6), or an 
error diffusion process. The dithering functions similar to multilevel halftoning 

25 algorithm and eliminates visual contours. 

The dot density and dot size controls are selected to obtain the desired gray 
level, e.g. tone curve, for each input value x(m, n). Advantageously, the dot size 
control and dot density control may be precomputed, and in one embodiment are 
designed with respect to each other and the characteristics of the printing system in a 

30 closed loop measurement process to optimize the quality of the image while matching 
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the printing device characteristics. Because the tone curve is controlled through the 
combination of the dot density control and dot size control there is an additional 
degree of freedom that may be used to optimize a variety of printing attributes 
including print quality and/or print stability. In other words, it is possible to achieve 
5 most desired tone curves by either varying the size of the dots or the density of the 
dots. Therefore, a particular choice of a dot size necessitates a specific dot density to 
achieve the desired tone curve. The object is then to select the dot size, at each 
desired output density, that produces the "best" printed output. The dot density LUT 
122 and dot size LUT 124 are designed by experimentally measuring the quality of 

10 the printed halftones at each combination of dot size and dot density, and then 

choosing the best combination of dot size and density for each desired output gray 
level. Because the dot density control and dot size control are designed by directly 
measuring the print quality, the resulting design accounts for attributes of both the 
printer and the halftoning system used. 

15 Fig. 10 is a flow chart 180 showing a method of optimizing the dot size and 

dot density look-up tables using a closed loop measurement process in accordance 
with an embodiment of the present invention. As shown in Fig. 10, a test page is 
printed for each fixed dot size over the range of dot densities (block 1 82). Let n be 
the dot density, 0 be the dot size, and g be the dot gray level. The dot gray level g is 

20 measured in units of gamma corrected absorbtance ranging from 0 to 1 . So the 

normalized luminance is given by / = (1 - g) y where y = 2.2. For each dot density and 
dot size, the output absorptance, i.e., gray level ge(n), and the print distortion, D e (n) 
are measured from the printed test pages that are scanned at 600 dpi (block 184). 
Each test page consists of an array of 16 x 16 blocks printed with a fixed dot size and 

25 dot densities ranging from 0 to 255 where 255 corresponds to all dots on. The blocks 
are of size 0.35in. x 0.35in and placed in random order to avoid systematic 
measurement errors when printing across a page. 

To measure the dot distortion D e (n) for each dot size 0 and dot density n, a 
frequency weighted total squared error is used as the metric, where the frequency 

30 weighting is chosen to approximate the human visual system (HVS) response. The 
HVS model is a linear shift-invariant low pass filter. The frequency response of this 
filter is given by 
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aL b expC-^+v 2 ) 

v) = v equ. 5 

clnL + rf 

where L is the average luminance in cd/m^ a = 131. 6, b = 0.3188, c = 0.525 and d = 
3.91. 

Let f(m, n) denote the measured gray leve/ measured in units of gamma 
5 corrected absorptance so that I(m,n) = (1- f(m, n)) y where I(m,n) is the normalized 
absorptance and y = 2.2. Then let 

£ = ^7 X f( m > n ) ec l u ' 6 

™ m,n 

be the average absorptance level. The windowed error is computed as: 

e(m, n) = (f (m, n) - g)w(m, n)(4) equ. 7 

10 where w(m, n) is a Hamming window. The discrete Fourier transform (DFT) of 
e(m,n) is calculated to form E(k, /), and the print distortion metric is then given by 



D = 



J^Xi^'Oi 2 H &>0 equ. 8 



where H(k, /) is the appropriately sampled version of H(u, v). 

Figs. 1 1 and 12 are graphs 200 and 210 showing the measured tone curve, 

15 ge(n), i-e. ? the output absorptance, and the print distortion, De(n), respectively, versus 
input dot density for a fixed value of dot size 9 = 1, where a square dot shape is used 
as discussed in reference to Fig. 9. As shown in Figs. 1 1 and 12, the measurements of 
the tone curve and the distortion are noisy, so they must be appropriately smoothed. 
Smoothing the tone curve also ensures that the smoothed tone curve is monotonically 

20 increasing. Smoothing is performed using a combination of median and linear 

filtering and polynomial curve fitting, which are well known in the art. For example, 
the raw data of the tone curve may be smoothed with a median filter of length 3 that is 
applied 5 times followed by a linear filter that is 3 times with a filter kernel [0.25, 0.5, 
0.25]. A 5 th order polynomial curve fitting is then applied to obtain a smooth tone 

25 curve. The printing distortion curve, for example, may be smoothed by squaring the 
measured data, applying a median filter of length 9, applying a linear filter 40 times 
with a filter kernel [0.25, 0.5, 0.25], and then taking the square root of the result. The 
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tone-corrected ramp for each dot size is printed out to see if the resulting curve is 
reasonable. Of course, any desired method of smoothing may be used to generate a 
reasonably smooth curve. 

Once the tone and distortion curves are determined for each dot size, e.g., 9 = 
5 0.3 to 1.4, these curves are used to determine the distortion at each gray level (block 
186). Each tone curve is inverted to compute the value of the dot density required to 
produce each gray level. 

n eis) = maxfe*(w) ^ g) equ. 9 

n 

Fig. 13 is a graph 220 showing a set of inverse tone curves for dot sizes 
10 ranging from 0.3 to 1.4. The x-axis is the desired output absorptance, while the y-axis 
is the dot density required to produce that absorptance level. 

The inverse tone curves are used to determine the print distortion as a function 
of gray level. 

D 0>g =D 0 (n e (g)) equ. 10 

15 Fig. 14 is a graph 230 showing a set of plots of the print distortion versus the 

desired output absorptance for dot sizes ranging from 0.3 to 1.4, where each curve 
corresponds to De )g as a function of g for fixed 9. These curves can be used to 
determine the best value of dot size for each desired output absorptance. For 
example, at the low values of absorptance, the smaller dot sizes generally produce 

20 lower distortion. 

With these measurements, the optimized dot size LUT 124 and dot density 
LUT 122 can be calculated (blocks 188 and 190). These two look-up tables denoted 
by 9j and n*, respectively, where i ranges from 0 to 255. The look-up table functions, 
gi and 0i may be optimized so that both are smooth and achieve nearly optimal print 

25 quality at each desired absorptance level, by minimizing a cost function shown below 
as a function of the look-up tables 9i. 

255 i 255 

cos/ = AZW-*") 2 equ. 11 

Here a is a parameter that controls smoothness of the result by forcing the values of 
the parameters Q\ to only change in small amounts. Because the print distortion 
30 function is only measured for discrete values of the dot size 9, the function must be 
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interpolated using e.g., cubic spline interpolation, for missing values. Coordinate 
decent optimization is used to minimize the cost function; however, other 
optimization methods can be used and may be desirable, particularly ones that are 
robust to local minima in the cost functional. 
5 Fig. 15 is a graph 240 showing a dot size curve for different a values as a 

function of desired absorptance, where the optimized dot size LUT 124 curve 242, 9i 
is a = 0.05. The smooth curve 242 is used as the dot size LUT 124 for the AM/FM 
halftoning system 120, shown in Fig. 5. Smaller values of a can be used to increase 
the smoothness of the dot size LUT, and larger values of a can be used to reduce the 
10 smoothness of the LUT. 

Once the look-up table Q[ is determined, the corresponding dot density LUT 
122 ni may be determined by using the inverse tone curves, n e (g) shown in graph 220 
in Fig. 13, so that: 

n t =n 0 (i) equ. 12 

15 Again, because n 6 (i) is not measured for each value of dot size 0, the intermediate 
values may be interpolated using, e.g., cubic spline interpolation. 

Fig. 16 is a graph 250 showing the dot density curve that corresponds to the 
curve 242 shown in graph 240 in Fig. 15 and is a function of the desired absorptance. 
The dot density curve shown in graph 250 is used as the dot density LUT 122 for the 
20 AM/FM halftoning system 120, shown in Fig. 5. 

While graphs 240 and 250 may be used as the dot size LUT 124 and dot 
density LUT 122 in the AM/FM halftoning system 120, in practice, the best results 
may require some manual adjustment of these curves. For example, the dot density 
curve in graph 250 is not monotone. While the resulting gray levels of the AM/FM 
25 halftoning system 120 are still monotone due to the increasing dot size at those gray 
levels, in practice, it may be more desirable to constrain the solution to have 
monotone increasing dot density. 

Figs. 17 and 18 are graphs 260 and 270, which show a dot size curve and dot 
density curve, respectively, as used with an eight-bit PWM code embodiment in 
30 which the pixel pairs are used, as described in reference to Figs. 7 A and IB. 

It should be understood that variations of the present AM/FM halftoning 
system are possible. For example, Fig. 19 shows an embodiment of the AM/FM 
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halftoning system 300, which is similar to the embodiments shown in Fig. 3, like 
designated elements being the same. AM/FM halftoning system 300, however, 
includes a dot size diffusion unit 302, which performs error diffusion for dot size in a 
manner similar to the dot density diffusion system described above, e.g., in reference 
5 to Fig. 6. The dot size diffusion unit 302 operates with dot size modulation unit 110 
within the AM modulation unit 304. This embodiment may be particularly useful 
when the PWM levels are limited. Thus, for example, if the PWM code is a two-bit 
code, four PWM levels are possible. With only four PWM levels, contouring artifacts 
or discontinuities will arise. With the use of dot size diffusion unit 302, the 

10 contouring artifacts may be minimized. 

The source code in Appendices Bl and B2 written in C programming 
language are an example of an implementation of the present invention using a two- 
bit PWM code to control a printer system having a pulse width modulated laser 
capability for each printed pixel. Appendix B3 lists the particular coefficients used in 

15 the Tone-Dependent Error Diffusion process for the two-bit code used in conjunction 
with Appendices Bl and B2. 

Fig. 20 shows a block diagram illustrating another possible embodiment of the 
present AM/FM halftoning system. Fig. 20 shows an integrated AM/FM halftoning 
system 350 in which the dot density control is integrated into the dispersed dot 

20 halftoning system 352 and the dot size control is integrated into the dot size 

modulation system 354. It should be understood, of course, that other variation exist, 
such as integrating only the dot density control into the dispersed dot halftoning 
system 352 or only the dot size control into the dot size modulation system 354. The 
integrated AM/FM halftoning system 350 has potential advantages of reducing 

25 quantization artifacts, reducing complexity, and improving halftone quality. 

In one embodiment, the look-up tables and method for controlling the error 
diffusion process are provided on a computer readable medium, such as a 
microdiskette or floppy diskette as a printer driver. This printer driver is then 
installed into the computer so that the program is installed in the computer's RAM. 

30 Such a program may be also installed in the printer and, in one embodiment, installed 
in firmware within the printer. All logic functions may be implemented in hardware 
or software. Thus, for example, the image forming system may include a computer 
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that is coupled to the printing device, where a computer program executed by the 
computer includes instructions for implementing the functions of the present 
invention. The computer may be, e.g., a host computer, microprocessor or any other 
appropriate device. If hardware is used, the various table values would be available to 
5 the circuitry implementing the halftone method via bus lines. The method may also 
be carried out by an ASIC, which controls the timing and transfer of data to the 
various logic devices and look-up tables as well as to and from the image map, as 
would be understood by those skilled in the art after reading this disclosure. 

While particular embodiments of the present invention have been shown and 
10 described, it will be obvious to those skilled in the art that changes and modifications 
may be made without departing from this invention in its broader aspects and, 
therefore, the appended claims are to encompass within their scope all such changes 
and modifications as fall within the true spirit and scope of this invention. 
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CLAIMS 



What is claimed is: 



5 



1. 



A method of halftoning an image, said method comprising: 
modulating the dot density of said image; and 
modulating the dot size of printed dots to obtain a printed halftone 



image. 



10 



2. 



The method of Claim 1, further comprising: 

inputting an input pixel value for a pixel location; 

wherein said modulating the dot density modulates the spacing of the 
dot to be printed for said pixel location with respect to preceding and 



subsequent pixel locations and said modulating the dot size modulates the size 
of said dot to be printed for said pixel location. 

3 . The method of Claim 1 , further comprising: 



controlling the dot density for said pixel location within said image using said 
input pixel value and performing dispersed dot halftoning to produce a dot 
position based on said dot density. 

4. The method of Claim 3, wherein said modulating the dot size of said image 
20 comprises controlling the dot size for said pixel location within an image using said 

input pixel value and performing dot size modulation based on said dot size and said 
dot position. 

5. The method of Claim 4, wherein controlling the dot density for said pixel 
location comprises generating a dot density value based on said input pixel value, said 

25 dot density value being used to perform said dispersed dot halftoning, and wherein 
said controlling the dot size for said pixel location comprises generating a dot size 
value based on said input pixel value, said dot size value being used to perform dot 
size modulation. 



15 



inputting an input pixel value for a pixel location within said image; 
wherein modulating the dot density of said image comprises 
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6. The method of Claim 5, wherein said generating a dot density value is 
performed using a first look up table and generating a dot size value is performed 
using a second look up table. 

7. The method of Claim 4, wherein controlling the dot density for said pixel 

5 location and performing dispersed dot halftoning are performed in a first integrated 
process and controlling the dot size for said pixel location and performing dot size 
modulation are also performed in a second integrated process. 

8. The method of Claim 4, wherein said modulating the dot size of said image 
further comprises diffusing at least a portion of the dot size for the pixel location to at 

10 least one subsequently processed pixel location. 

9. The method of Claim 1, wherein said modulating the dot density is performed 
using at least one of error diffusion, dispersed dot screening, and iterative search 
based halftoning. 

10. The method of Claim 2, wherein said modulating the dot density is performed 
15 using tone dependent error diffusion. 

1 1 . The method of Claim 10, further comprising: 

generating a dot density value based on said input pixel value, said dot 
density value being used in said tone dependent error diffusion; 
said tone dependent error diffusion, comprising: 
20 combining said dot density value with at least one previous 

error value to produce a modified pixel value; 

comparing said modified pixel value with a threshold value to 
produce a halftone value for said pixel location; and 

using said halftone value for said pixel location to produce an 
25 error value that is diffused to at least one subsequently processed pixel. 

12. The method of Claim 1 , wherein said modulating the dot size is performed 
using pulse width modulation. 
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13. The method of Claim 12, wherein said pulse width modulation provides a first 
value indicating the width of the pulse and a second value indicating the justification 
of said pulse within a pixel location. 

14. The method of Claim 13, wherein when a printed dot is larger than one pixel, 
5 the pulse within adjacent pixels are justified together so that the dot may be formed 

with one continuous pulse. 

A method of printing a halftone image, said method comprising: 

inputting an input pixel value for a pixel location within an image to be 
printed; 

generating a dot density control value using said input pixel value for 
said pixel location; 

producing a halftone value for said pixel location using said dot 
density control value; 

generating a dot size control value using said input pixel value for said 
pixel location; and 

producing a size modulated halftone value for said pixel location based 
on said dot size control value and said halftone value. 

16. The method of Claim 15, wherein said producing a halftone value for said 
pixel location also uses accumulated errors diffused from at least one other pixel 

20 location and provides an error for said pixel location to be diffused to at least one 
subsequently processed pixel location. 

17. The method of Claim 15, wherein said producing a halftone value for said 
pixel location is performed using tone dependent error diffusion. 

18. The method of Claim 15, wherein generating a dot density control value is 
25 performed using a first look up table and generating a dot size control value is 

performed using a second look up table. 

19. The method of Claim 15, wherein said producing a size modulated halftone 
value is performed using pulse width modulation. 
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20. The method of Claim 19, wherein said pulse width modulation is performed 
using a look up table. 

21 . The method of Claim 1 5, further comprising diffusing at least a portion of said 
dot size value for the pixel location to at least one subsequently processed pixel 

5 location. 

22. A method of optimizing a dot size look-up table and a dot density look-up 
table for a printing system that uses dot size modulation and dot density modulation, 
the method comprising: 

printing at least one test page showing the combinations of dot sizes 
10 and dot densities; 

measuring the output absorptance for each combination and the print 
distortion for each combination; 

determining the print distortion at each output absorptance using the 
measured output absorptance and the measured print distortion for each 
15 combination; 

calculating the optimized dot size look-up table using the print 
distortion at each output absorptance; and 

calculating the optimized dot density look-up table using the print 
distortion at each output absorptance. 

20 23 . The method of Claim 22, wherein determining the print distortion comprises: 
inverting the output absorptance for each combination to compute the 
value of the dot density required to produce each output absorptance; and 

using the inverted output absorptance and the print distortion for each 
combination to determine the print distortion as a function of output 
25 absorptance. 

24. The method of Claim 22, wherein calculating the optimized dot size look-up 
table comprises minimizing a cost function of the print distortion at each output 
absorptance as a function of dot size. 
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25. The method of Claim 24, wherein calculating the optimized dot density look- 
up table further comprises using the calculated optimized dot size look-up table. 

26. A printing system including a printer and a computer, the printing system 
printing a dot density and dot size modulated image, the method performed by the 

5 printing system comprising: 

performing dispersed dot halftoning for a pixel location based on an 
input pixel value for said pixel location; 

performing dot size modulation for said pixel location based on the 
results of said dispersed dot halftoning and said input pixel value. 

10 27. The image forming device of Claim 26, wherein the printer is one of an 
electrophotographic printing device, electrophotographic copying device, and an 
inkjet printer. 

28. An image forming system comprising: 

a computer; 

15 a printing device coupled to said computer; 

a computer program executed by said computer, wherein said 
computer program comprises computer instructions for: 
receiving an image; 

modulating the dot density of said image; and 
20 modulating the dot size of printed dots to obtain a 

halftone for said image. 

29. The image forming system of Claim 28, wherein: 

receiving an image comprises receiving an input pixel value for a pixel 
location within said image; 
25 modulating the dot density of said image comprises generating a dot 

density control value using said input pixel value for said pixel location and 
producing a halftone value for said pixel location using said dot density 
control value; and 

modulating the dot size comprises generating a dot size control value 
30 using said input pixel value for said pixel location and producing a size 
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modulated halftone value for said pixel location based on said dot size control 
and said halftone value. 

30. The image forming system of Claim 29, wherein said producing a halftone 
value for said pixel location comprises adding at least a portion of accumulated errors 

5 from at least one other pixel location with said dot density control value and diffusing 
an error for said pixel location to at least one subsequently processed pixel location. 

3 1 . The image forming system of Claim 28, wherein said computer is a 
microprocessor. 

32. The image forming system of Claim 28, wherein said image forming system is 
10 one of an electrophotographic printing device, electrophotographic copying device, 

and an inkjet printer. 

33. An image forming device comprising: 

a dot density control unit; 
a dot size control unit; 

15 a dispersed dot halftoning unit coupled to said dot density control unit; 

and 

a dot size modulation unit coupled to said dispersed dot halfoning unit 
and said dot size control unit; and 

a printing device coupled to said dot size modulation unit. 

The image forming device of Claim 33, wherein: 

said dot density control unit is configured to receive an image and 
produce a dot density value for at least a portion of said image; 

said dot size control unit is configured to receive said image and 
produce a dot size value for said at least a portion of said image; 

said dispersed dot halftoning unit is configured to receive said dot 
density value and produce a halftone value for said at least a portion of said 
image; 

said dot size modulation unit is configured to receive said halftone 
value and said dot size value and produces an output signal representing said 



20 34. 



25 
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halftone value and a modulated size value for said at least a portion of said 
image; 

said printing device is configured to receive said output signal and 
print said at least a portion of said image. 

5 35. The image forming device of Claim 33, wherein said dispersed dot halftoning 
unit is a tone dependent error diffusion unit. 

36. The image forming device of Claim 33, further comprising a dot size diffusion 
unit coupled to said dot size modulation unit. 

37. The image forming device of Claim 33, wherein said dot density control unit 
10 is integrated into said dispersed dot halftoning unit. 

38. The image forming device of Claim 33, wherein said dot size control unit is 
integrated into said dot size modulation unit. 
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ABSTRACT 

A combined dot density and size modulation system uses dispersed dot 
halftoning in conjunction with dot size modulation to produce a halftone image in 
which both the density and size of the dots are modulated to control overall gray level. 
5 The dot density and size modulation system offers advantages over pure dot density 
modulation systems or pure dot size modulation systems because it allows an extra 
degree of flexibility which can be used to increase the visual quality of the halftoned 
pattern and/or increase the robustness of the halftoning to printer artifacts and 
variations. An input pixel value is used to independently produce a dot density value 

10 and a dot size value. The dot density value and dot size values may be obtained from, 
e.g., look up tables that have been optimized for print quality and printer stability. 
Dispersed dot halftoning is used to provide a halftone value for the desired pixel 
location using the dot density value. The dispersed dot halftoning may be, e.g., tone 
dependent error diffusion. The halftone value and the dot size value for the pixel 

15 location is then used to generate a modulated code, e.g., a pulse width modulated 
code, to the printer. The modulated code may include both the pulse width of the 
desired dot for the pixel location as well as the justification, e.g. left, center, or right, 
for the pixel location. The dot density and size modulation system is particularly 
useful in modern electrophotographic printing systems that allow the printed dot size 

20 to be almost continuously varied through the specification of a pulse width 
modulation (PWM) code. 
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COMBINED DOT DENSITY AND DOT SIZE MODULATION 
Zhen He 

Charles A Bouman 
Qian Lin 

10003284 
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/* amfm_8bit.c file */ 

/* 8 bits/pixel am/fm halftoning algorithm (with partial doting) */ 
/* Image length and width are assumed to be multiple of 8 */ 
/* One row serpentine TDED with suppressing each other dot */ 
/* The output is a tiff file containing 8bit pwm codes */ 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <time.h> 
#include " coef . h " 
#include "tiff .h" 
#include "allocate . h" 

void amfm (unsigned int, unsigned int, unsigned char *, unsigned char *,\ 

short *, unsigned char * * , TDED PARA *, short * , short *); 
void amfm_ed_8bits (unsigned int, unsigned int, unsigned char **, unsigned char**,\ 

unsigned char * * , TDEDPARA * f short *, short *); 
int main (int argc, char ** argv) 

{ 

int i , j ; 

unsigned int height , width; 
FILE * fp; 

struct TIFF_img input_img, output_img, mid; 
time_t first, second; 

TDEDPARA *tdedpara = &TDEDcoef f [0] ; 

short *dotdensityLUT = &OptDensityLUT [0 ] ; 

short *dotsizeLUT = &OptSizeLUT [0 ] ; 

if(argc<3) { 

printf ( "usage : %s input_img.tif output_img . tif \n" , argv [0] ) ; 
return 1 ; 

} 

/* read the input image */ 

if ( (fp=f open (argv [1] , " rb " ) ) ==NULL) { 



printf {"can not open file %s l\n" , argv [1] ) ; 
exit (2) ; 

} 

if (read__TIFF (fp, &input_img) ) { 

printf ( "error reading input file\n"); 
exit (3) ; 

} 

f close (f p) ; 

if ( (fp=fopen("dbshalf . tif " , "rb") ) ==NULL) { 

fprintf (stderr, "can not open file: dbshalf . tif \n" ) ; 
exit (1) ; 

} 

if (read__TIFF (fp, &mid) ) 
{ 

fprintf (stderr , "error reading file\n" ); 
exit (1) ; 

} 

f close (f p) ; 

/* Set variable to do timing of algorithm */ 
first = time (NULL) ,- 

/* Modify image width to make sure each strip is multiple of 8 */ 
width = f loor ( input_img . width/ 8) *8 ; 
height = f loor (input_img- height/ 8) *8 ; 

/* Allocate memory for entire fm output image. */ 
get_TIFF( &output_img , height, width, ' g* ) ; 

amfm_ed__8bits (height , width, input_img .mono, out put _img .mono, \ 
mid . mono , tdedpara , dot dens ityLUT, dotsizeLUT) ; 

/* show the run time */ 
second = time (NULL) ; 

fprintf (stdout, "\nFinished AM/ FM and writing results . \n" ) ; 
fprintf (stdout, "Cum. run time: %f sec . \n" , diff time (second, first ) ) 

/* write PWM codes image */ 

if( (fp = f open (argv [2] , "wb") ) ==NULL) { 

printf ("cannot open file %s\n'\ argv[2]); 

exit (4) ; 

} 

if (write_TIFF (fp, &output__img) ) { 

printf ("\nError writing TIFF file %s\n" , argv[2]); 
return 1 ; 

} 

f close (fp) ; 

/* free the space */ 

f ree_TIFF (& (output_img) ) ; 

f ree_TIFF (& (input_img) ) ; 

f ree_TIFF (&(mid) ) ; 

f flush (stdout) ; 

return 0 ; 

} 



void amf m_ed__8bits ( 

unsigned int height, /* 
unsigned int width, /* 
unsigned char ** contone_img , 



unsigned char ** token_img, 
unsigned char ** dbs_screen 
part */ 

TDEDPARA *tdedpara, 
short *dotdensityLUT, 
short *dotsizeLUT) 



{ 



I* 
/* 
/* 



Input image height */ 
Input image width */ 

/* Input image [height] [width] */ 
Output token image [height] [width] */ 

/* DBS screen used in thresholding of fm 

Tone -dependent error diffusion parameters */ 
Optimal dot density curve */ 
Optimal dot size curve */ 



short *fm__err; 
unsigned int 

/* initialize first row of fm error buffer */ 
srand(l) ; /* fix the seed */ 

fm_err = (short*) malloc (sizeof (short ) * (width+2) ) ; 
for(j = 0; j<width+2; 

fm_err[j] = (rand ( ) %128 -64 ) ; /* initialization */ 

/* Process the input image with 2 rows each time */ 
for(i=0; i<height; i+=2) { 

if((i%600) == 0) print f ( "amf m_ed: starting row %d\n" , i) ; 

amfm (width, i, contone_img [i] , token_img [i] , f m_err , dbs__screen, \ 
tdedpara, dotdensityLUT, dotsizeLUT) ; 

} 

free (fm_err) ; 
return; 



> 

/* This subroutine only processes 2 rows */ 
/* Assume width of image is multiple of 8 */ 
void amfm ( 

unsigned int width, /* Input image width */ 
unsigned int i, /* ith row */ 

unsigned char *img__in, /* ith row of input image array */ 

unsigned char *img_out, /* ith row of output image array */ 

short *fm_err, /* FM error buffer */ 

unsigned char ** dbs_screen, /* dbs_screen [SCREENHEIGHT] [SCREENWIDTH] 

TDEDPARA *tdedpara, /* Tone -dependent error diffusion parameters */ 
short * dot dens ityLUT, /* Optimal dot density curve */ 
short *dotsizeLUT) /* Optimal dot size curve */ 



{ 

short fm__tmp , thresholding ; 
short *fm_err _j?tr , *tded_ptr ; 
short pixela, pixelb, output; 
int j ; 

unsigned char *img__in_j?tr , * img_out__pt r , *dbs_pat_rowptr ; 

short dotdensity, mod_input , error; 

short Wl, W2, W3, W4 , T2 , DT , el, e2 , e3, e4 ; 

FILE *fp; 

/* */ 



/* serpentine even rows */ 
/* */ 

/* initial points */ 
f m__tmp = 0 ; 

fm_err_ptr = fm_err+l; 
img_in_ptr = img_in; 
img_out_jptr = img_out ; 

dbs__jpat_rowptr = dbs_screen [ (i++) %SCREENHEIGHT] ; /* Inverse dbs pattern */ 

/* Index through pixels in pairs */ 
for(j = 0; j<width; j = j +2 ) { 

/* First process FM (dot density) for left pixel in pixel pair. */ 

/* Get first pixel */ 
pixela = * (img_in_jptr++) ; 

/* Use look-up-table to get dot density */ 
dotdensity = dotdensityLUT [pixela] ; 

/* Compute look-up table entries for tone dependent error diffusion */ 
tded_jptr = (short*) (tdedpara + dotdensity) ; 
T2 = * (tded_ptr++) 
DT = * (tded_jptr++) 
Wl = * (tded_jptr++) 
W2 = * (tded_jptr++) 
W3 = * (tded_jptr++) 
W4 = *tdedj>tr ; 

/* compute dotdensity modified by diffused error */ 
mod_input = dotdensity + *fm_errj)tr; 

/* Threshold modifed dotdensity */ 

thresholding = mod_input - (dbs_pat_rowptr [ j %SCREENWIDTH] * DT + T2) ; 
output = (thresholding > 0) ? 255 : 0; 

/* Compute weighted errors */ 
error = output - mod_input; 
el = (Wl * error) »8; 
e2 = (W2 * error) >>8; 
e3 = (W3 * error) >>8; 
/*e4 = (W4 * error) >>8;*/ 
e4 = error - el - e2 -e3; 

/* Diffuse error forward in 1-D error buffer */ 
* ( - -fm_errj)tr) -= e4 ; 
* (++fm_err_ptr) = fm_tmp - e3 ; 
* (++f m_err_ptr) -= el; 
f m_tmp - -e2; 

/* Now process FM (dot density) for right pixel in pixel pair. */ 
/* Use same TDED parameters as for left pixel. */ 

/* Get second pixel */ 
pixelb = * (img_in_j)tr+ + ) ; 



/* Use look-up-table to get dot density */ 



dotdensity = dotdensityLUT [pixelb] ; 
mod_input = dotdensity + *fm_errjptr; 

error = - mod_input; /* suppress dot firing at this pixel */ 

el = (Wl * error) >>8; 
e2 = (W2 * error) >>8; 
e3 = (W3 * error) >>8; 
/*e4 = {W4 * error) >>8;*/ 
e4 = error - el - e2 -e3; 
/* Using the tded weights 
* ( - - f m_err_jptr) - - e4 ; 
*(++fm_err_ptr) = fm_tmp 
* (++fm_err_ptr) - = el; 
fm_tmp = -e2; 

* Begin section on dot size rendering with partial doting */ 

if (output) { 

/* Left pixel */ 

* (img__out_ptr+ + ) = (dotsizeLUT [pixela] >>1) +NE WRIGHT ; 
/* Right pixel */ 

if (dotsizeLUT [pixela] & 1) /* Take care of quantization error */ 

* (img_out_ptr++) - ( (dotsizeLUT [pixelb] +1) >>1) + NEWLEFT ; 
else 

*(img_outjptr++) = (dotsizeLUT [pixelb] > >1 ) + NEWLEFT; 

} 

else { 

* (img_out_ptr++) = NEWRIGHT; 
* (img_out_ptr++) = NEWLEFT ; 

} 

} /* end of ith row */ 

* */ 

* serpentine odd rows */ 
* */ 

f m_tmp = 0 ; 

/* Set fm error buffer pointer to the end of fm_err buffer */ 
fm_err__ptr = fm_err+ width - 1; /* offset by 1 */ 
img_in_ptr = img_in+width*2 -2 ; 
img_out_ptr = img_out+width*2-l ; 

*img_out_ptr = NE WRIGHT ; /* Take care of the last pixel in odd row */ 
img_out_ptr -= 2; 

dbs j»at_rowptr = dbs_screen [i%SCREENHEIGHT] ; 

/* Index through pixels in pairs */ 
for(j = width-2; j> 0; j=j-2) { 

/* First process FM (dot density) for right pixel in pixel pair */ 

/* Get right pixel */ 

pixela = * (img_in__ptr--) ; 

/* Use look-up- table to get dot density */ 
dotdensity = dotdensityLUT [pixela] ; 



of the left pixel */ 
- e3; 



/* Compute look-up table entries for tone dependent error diffusion 



tded_j>tr = (short*) (tdedpara + dotdensity) ; 



T2 


= * (tded_ptr++) ; 


DT 


= * (tded_ptr++) ; 


Wl 


= * (tded_jptr++) ; 


W2 


= * (tded_ptr++) ; 


W3 


= * (tded_ptr++) ; 


W4 


= *tded_ptr; 



/* Compute dotdensity modified by diffused error */ 
mod_input = dotdensity + *f m_err_j?tr ; 

/* suppress this dot and compute the error */ 
error - - mod_input ; 

/* Compute weighted errors */ 

el = (Wl * error) >>8; 

e2 = (W2 * error) >>8; 

e3 = <W3 * error) >>8; 

/*e4 = ( {W4 * error) >>8) ; */ 

e4 = error - el - e2 -e3; 

/* duffuse error forward in 1-D error buffer */ 
* (++fm_err_ptr) -= e4 ; 

* (--fm_err_ptr) = f m_tmp - e3 ; 
* ( - -f m_err_ptr ) - = el; 
fm_tmp = -e2; 

/* Now process FM (dot density) for Left pixel in a pair */ 

/* Get second pixel */ 
pixelb = * ( img__in_ptr- - ) ; 

/* Use look-up-table to get dot density */ 
dotdensity = dotdensityLUT [pixelb] ; 

mod_input = dotdensity + *fm_errj>tr; 

/* Threshold modi fed dotdensity */ 

thresholding - mod_input - (dbs_jpat_rowptr [ ( j -1) %SCREENWIDTH] * DT + T2 ) 
output = (thresholding > 0) ? 255 : 0; 

error = output - mod_input ; 

el = (Wl * error) >>8; 
e2 = (W2 * error) >>8; 
e3 = (W3 * error) >>8; 
/*e4 = <W4 * error) >>8;*/ 
e4 = error - el - e2 -e3; 

* (++fm_err_ptr) -= e4 ; 

* { - -f m_errj>tr) = f m_tmp - e3 ; 
* ( - -f m_err_ptr ) - = el; 
fm__tmp = -e2; 

/* Begin section on dot size rendering with partial doting */ 
if (output) { 



/* Left pixel */ 

* (img_out_ptr++) - (dotsizeLUT [pixelb] »1) + NEWRIGHT; 
/* Right pixel */ 

if (dotsizeLUT [pixelb] & 1) /* Take care of quantization error */ 

*img_out_jptr = ( {dotsizeLUT [pixela] +1) »1) + NEWLEFT; 
else 

*img_outj>tr - {dotsizeLUT [pixela] »1) + NEWLEFT; 

} 

else { 

* (img_out_ptr++) = NEWRIGHT; 
*img_out_ptr = NEWLEFT ; 

} 

img__out_ptr - = 3; 

* (++img_out_ptr) = NEWLEFT; /* Take care of the first column */ 
return; 

} 
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/* coef.h file */ 



#define SCREENHEIGHT 128 
#define SCREENWIDTH 12 8 



#define NE WRIGHT OxcO 
#define NEWLEFT 0x4 0 
#define NEWCENTER 0x00 



#define Fl 0x0007 
#define F2 0x0003 
#define F3 0x0005 
#define F4 0x0001 



/* Floyd-Steinberg 

/* Floyd-Steinberg 

/* Floyd-Steinberg 

/* Floyd-Steinberg 



Weights 7/16 in Q4 

Weights 3/16 in Q4 

Weights 5/16 in Q4 

Weights 7/16 in Q4 



typedef struct TDEDPARA 

{ 

short T2; 
short DT; 
short Wl; 
short W2; 
short W3; 
short W4; 
} TDEDPARA; 



stat 


ic 


TDEDPARA 


TDEDcoef f [129] 


{76, 


0, 


181, 


0, 


3, 72}, 


{76, 


0, 


181, 


o, 


3, 72}, 


{79 , 


0, 


172 , 


1, 


2, 81}, 


{80, 


o, 


161, 


14, 


18, 63}, 


{82, 


0, 


159, 


1/ 


37, 59}, 


{83 , 


0 , 


149, 


6, 


5, 96}, 


{83, 


0, 


141, 


30, 


0, 85}, 


{85, 


o, 


138, 


13 , 


0, 105}, 


{86, 




144 , 


10, 


1, 101}, 


{85, 


0, 


129, 


48, 


3, 76}, 


{86 , 


o, 


123 , 


31, 


1, 101}, 


{87 , 


o, 


123, 


29, 


3, 101}, 



={ 



87, 0, 115, 28, 5, 108} 
89, 0, 138, 19, 18, 81} 
89, 0, 111, 17, 51, 77} 

88, 0, 115, 31, 0, 110} 

87, 0, 120, 16, 16, 104 

88, 0, 139, 12, 0, 105} 

89, 0, 122, 19, 17, 98} 

90, 0, 112, 32, 0, 112} 

91, 0, 98, 34, 20, 104} 

90, 10, 123, 16, 26, 91 
93, 8, 126, 1, 74, 55}, 
89, 10, 89, 26, 71, 70} 
89, 10, 89, 22, 43, 102 
89, 12, 91, 21, 34, 110 
88, 12, 85, 24, 30, 117 
88, 14, 85, 23, 30, 118 
84, 24, 113, 27, 13, 103 

82, 26, 113, 33, 0, 110 

83, 26, 109, 29, 9, 109 

84, 28, 106, 21, 29, 100 

85, 28, 103, 13, 56, 84 
96, 2, 102, 16, 57, 81} 

93, 6, 102, 25, 28, 101 

91, 12, 102, 24, 32, 98 

96, 2, 103, 24, 23, 106 

94, 10, 99, 17, 62, 78} 

95, 6, 110, 12, 110, 24 

97, 4, 114, 12, 112, 18 
97, 6, 114, 11, 113, 18 

96, 8, 111, 14, 110, 21 
94, 12, 102, 17, 109, 28 

94, 8, 79, 32, 108, 37} 

95, 6, 74, 35, 110, 37} 

97, 2, 70, 35, 111, 40} 
97, 4, 68, 33, 112, 43} 

97, 6, 69, 28, 112, 47} 

98, 6, 70, 22, 114, 50} 
97, 6, 68, 43, 113, 32} 
100, 4, 68, 22, 114, 52 

99, 6, 71, 24, 112, 49} 
102, 2, 70, 23, 114, 49 

100, 6, 68, 23, 114, 51 
100, 8, 66, 22, 116, 52 
100, 8, 66, 24, 116, 50 

96, 16, 75, 0, 122, 59} 
95, 16, 63, 0, 127, 66} 
95, 16, 56, 0, 130, 70} 

97, 14, 56, 0, 132, 68} 
97, 16, 59, 0, 132, 65} 

97, 16, 60, 0, 133, 63} 

98, 16, 62, 0, 133, 61} 
95, 26, 98, 0, 109, 49} 

97, 20, 65, 0, 132, 59} 

98, 18, 61, 0, 132, 63} 

99, 18, 63, 0, 131, 62} 

100, 16, 58, 0, 133, 65 
100, 16, 58, 0, 131, 67 



}. 

r 



{101, 16, 60, 0, 131, 65}, 
{101, 16, 63, 0, 129, 64}, 
{101, 16, 58, 0, 129, 69}, 
{102, 16, 71, 0, 123, 62}, 
{103, 8, 68, 23, 114, 5l}, 
{103, 8, 66, 22, 116, 52}, 
{105, 6, 68, 22, 115, 5l}, 
{106, 4, 70, 22, 114, 50}, 
{108, 2, 69, 23, 113, 5l}, 
{105, 8, 68, 22, 114, 52}, 
{108, 6, 70, 20, 115, 51}, 
{106, 8, 69, 27, 112, 48}, 
{109, 2, 65, 35, 112, 44}, 
{110, 4, 69, 34, 111, 42}, 
{110, 6, 72, 35, 110, 39}, 
{114, 0, 73, 34, 111, 38}, 
{110, 12, 94, 21, 108, 33}, 
{111, 12, 102, 15, 110, 29}, 
{116, 6, 114, 10, 113, 19}, 
{96, 16, 92, 16, 67, 8l}, 
{100, 12, 95, 17, 67, 77}, 
{101, 12, 97, 19, 67, 73}, 
{99, 4, 101, 20, 45, 90}, 
{93, 4, 103, 25, 25, 103}, 
{94, 8, 101, 25, 33, 97}, 
{78, 24, 99, 26, 19, 112}, 
{81, 26, 104, 22, 24, 106}, 
{82, 26, 102, 26, 25, 103}, 
{91, 26, 109, 14, 46, 87}, 
{104, 10, 82, 0, 95, 79}, 
{107, 8, 83, 0, 97, 76}, 
{105, 8, 87, 2, 84, 83}, 
{81, 14, 86, 27, 25, 118}, 
{99, 12, 122, 0, 37, 97}, 
{102, 10, 117, 0, 45, 94}, 
{103, 10, 90, 21, 64, 81}, 
{105, 12, 122, 4, 51, 79}, 
{101, 12, 126, 9, 29, 92}, 
{88, 12, 121, 25, 0, 110}, 
{85, 12, 114, 25, 1, 116}, 
{89, 10, 109, 23, 10, 114}, 
{86, 12, 112, 29, 1, 114}, 
{89, 12, 119, 31, 0, 106}, 
{94, 10, 123, 37, 1, 95}, 
{93, 8, 117, 63, 1, 75}, 
{99, 6, 118, 75, 9, 54}, 
{97, 6, 120, 43, 3, 90}, 
{ill, 6, 121, 35, 32, 68}, 
{95, 6, 116, 54, 0, 86}, 
{107, 6, 125, 39, 15, 77}, 
{93, 34, 137, 27, 19, 73}, 
{85, 44, 139, 33, 16, 68}, 
{87, 48, 146, 31, 23, 56}, 
{87, 44, 148, 22, 10, 76}, 
{93, 40, 152, 22, 11, 7l}, 
{97, 44, 159, 4, 28, 65}, 
{95, 42, 161, 25, 4, 66}, 



{103, 48, 176, 3, 44, 33}, 
{101, 56, 165, 27, 55, 9}, 
{97, 56, 165, 27, 55, 9}, 

}; 
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/* amfm_2bit.c file */ 

/* 2 bits/pixel am/fm halftoning algorithm: part A */ 
/* Process input image in 4 strips */ 
/* Assume width of each stripe is multiple of 8 */ 
/* One row serpentine TDED */ 

/* One path amfm with dot size error diffusion */ 

/* A tiff image file containing bit-packed tokens is generated for amfm _jpwm. 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <time.h> 
#include "coef.h" 
#include "tiff.h" 
#include "allocate .h" 

void get_stripe {int r int, unsigned char **, int , int , unsigned char **) ; 
void cut_out_result (int , int , int, int , unsigned char **, unsigned char **) ; 
void amfm (unsigned int, unsigned int, unsigned char *, unsigned char * , short *, 

short *,unsigned char * * , TDED PARA * , TOKENLUT * , short * , short *) ; 
void amfm_ed_2bits (unsigned int, unsigned int,unsigned char **, unsigned char 

unsigned char **, TDED PARA *, TOKENLUT *, short *, short *); 
int main (int argc, char ** argv) 

{ 

int i, j , width, start_point, cut_off set, cut_width, store_of f set , stripe__width; 
FILE * fp; 

struct TIFF_img input_img, output_img, mid; 
time__t first, second; 

unsigned char **stripe, **output_stripe ; 
TDED PARA *tdedpara = kTDEDcoef f [0] ; 
TOKENLUT *tokenLUT = &TokenLUT [ 0 ] + 30; 
short *dotdensityLUT = &OptDensityLUT [0] ; 
short *dotsizeLUT = ScOptSizeLUT [0] ; 



if (argc<3 ) { 



printf ("usage: %s input_img.tif output_img . tif \n" , argv [0 ] ) ; 
return 1 ; 

} 

/* read the input image */ 

if ( (fp=f open (argv [1] , " rb" ) ) ==NULL) { 

printf ("can not open file %s l\n" , argv [1] ) ; 

exit (2) ; 

if (read_TIFF (fp, &input_img) ) { 

printf ( "error reading input file\n") ; 
exit (3 ) ; 

} 

f close (fp) ; 

if ( (fp=fopen( "dbshalf.tif", M rb") ) ==NULL) { 

f printf (stderr, "can not open file: dbshalf . tif \n" ) ; 
exit (1) ; 

} 

if (readJTIFF (fp, &mid) ) 
{ 

f printf (stderr , "error reading file\n"); 
exit (1) ; 

} 

f close (fp) ; 

/* Set variable to do timing of algorithm */ 
first = time (NULL) ; 

/* Modify image width to make sure each strip is multiple of 8 */ 
width = f loor (input_img .width/ 32 . 0) *32 ; 

/* Allocate memory for entire fm output image. */ 
get_TIFF( &output_img, input_img . height , width/4, 'g f ); 

/* Process 4 stripes independently */ 

cut_width = width/ 16; 

/* first stripe */ 

printf ( "\nProcess first stripe\n" ); 

stripe_width = width/4+OVERLAP_WIDTH/2 ; 

stripe = ( unsigned char ** 
)multialloc (sizeof (char) , 2 , input__img . height , stripe_width) ; 

output_s tripe = ( unsigned char ** 
Jmultialloc (sizeof (char) , 2 , input_img . height , stripe_width/4) ; 

start__point = 0; 

get_stripe (input_img . height , input_img . width, input_img . mono , start __point , stripe_wi 
dth, stripe) ; 

amfm_ed__2bits (input_img . height , stripe_width, stripe, output_stripe , \ 

mid.mono; tdedpara , tokenLXJT, dotdensityLUT, dotsizeLUT) ; 
cut_offset = 0; 
store__of f set = 0; 

cut_out__result ( input_img . height , cut_of f set , cut_width, store_offset, \ 

output_s tripe , output_img . mono) ; 
multifree (stripe, 2) ; 



multif ree (output^stripe , 2) ; 
printf ("\n") ; 

/* Second stripe */ 

printf ( "Process second stripe\n H ); 

stripe_width = width/4+OVERLAP_WIDTH; 

s tripe = { unsigned char ** 
) mult ialloc (sizeof (char) , 2 , input_img . height , stripe_width) ; 

output_s tripe = ( unsigned char ** 
)multialloc (sizeof (char) , 2 , input_img . height , stripe_width/4 ) ; 

Start_point = width/4 -0VERLAP_WIDTH/2 ; 

get_s tripe (input_img. height , input_img .width, input_img. mono, start_point , st ripe_wi 
dth, stripe) ; 

amfm_ed_2bits (input_img. height , stripe_width, stripe, output_stripe , \ 

mid. mono, tdedpara, tokenLUT, dot dens ityLUT, dotsizeLUT) ; 
cut_offset = 0VERLAP_WIDTH/8; 
store_offset = width/ 16; 

cut_out_result (input_img .height , cut_of f set , cut_width, store_of f set , \ 

output_s tripe , output_img . mono ) ; 
multif ree (stripe, 2) ; 
multif ree (output_stripe, 2) ; 
printf ("\n") ; 

/* Third stripe */ 

printf ("Process third stripe\n" ); 

stripe_width = width/4+OVERLAP_WIDTH; 

stripe = ( unsigned char ** 
) multialloc (sizeof (char) , 2 , input_img . height , stripe_width) ; 

output_stripe = ( unsigned char ** 
) multialloc (sizeof (char) , 2 , input_img . height , stripe_width/4) ; 

start_point = width/2 -OVERLAP_WIDTH/ 2 ; 

get_stripe (input_img . height , input_img . width, input_img .mono, start _jpoint , stripe_wi 
dth, stripe) ; 

amfm_ed_2bits (input_img . height , stripe_width, stripe, output_s tripe , \ 

mid. mono, tdedpara, tokenLUT, dot densityLUT, dotsizeLUT) ; 
Cut_offset = 0VERLAP_WIDTH/8; 
store_offset = width/ 8 ; 

cut_out_result (input_img . height , cut_of f set , cut_width, store_of f set , \ 

output__s tripe , output_img . mono) ; 
multif ree (stripe, 2) ; 
multifree(output_stripe,2) ; 
printf ("\n") ; 

/* Fourth stripe */ 

printf ( "Process fourth stripe\n"); 

stripe_width = width/4+OVERLAP_WIDTH/2 ; 

stripe = ( unsigned char ** 
) multialloc (sizeof (char) , 2 , input_img . height , stripe_width) ; 

output__s tripe = ( unsigned char ** 
) multialloc (sizeof (char) , 2 , input_img . height , stripe_width/4 ) ; 

start_point = width/4*3 -0VERLAP_WIDTH/2 ,- 

get_stripe (input_img . height , input_img . width , input_img .mono, start _point , stripe_wi 
dth, stripe) ; 

amfm_ed_2bits ( input_img . height , stripe_width , stripe , output__stripe , \ 



mid. mono, tdedpara, tokenLUT, dotdensityLUT, dotsizeLUT) ; 
cut_offset = 0VERLAP_WIDTH/8; 
store_offset = width/l6*3; 

cut_out_result ( input_img . height , cut_of f set , cut_width, store_of f set , \ 

output_s tripe , output__img . mono ) ; 
multifree (stripe , 2 ) ; 
multif ree (output_stripe , 2) ; 
printf ( f, \n n ) ; 



/* show the run time */ 
second = time (NULL) ; 

f printf (stdout, "\nFinished AM/FM and writing results . \n" ) ; 

f printf (stdout, "Cum. run time: %f sec . \n" , diff time (second, first) ) ; 

/* write PWM codes image */ 

if( (fp - fopen(argv[2] , "wb") )==NULL) { 

printf ("cannot open file %s\n" , argv[2]); 

exit (4) ; 

} 

if (write_TIFF (fp, &output__img ) ) { 

printf ("\nError writing TIFF file %s\n" / argv[2]); 
return 1 ; 

} 

f close (fp) ; 



/* free the space */ 

f reeJTIFF {& (output_img) ) ; 

free_TIFF (& (input_img) ) ; 

free_TIFF(&(mid) ) ; 

f flush (stdout) ; 

return 0; 



void amf m_ed_2bits ( 

unsigned int height, /* Input image height */ 

unsigned int width, /* Input image width */ 

unsigned char ** contone_img, /* Input image [height] [width] */ 

unsigned char ** token_img, /* Output token image [height] [width/4] */ 

unsigned char ** dbs__screen, /* DBS screen used in thresholding of fm 

part */ 

TDEDPARA *tdedpara, /* Tone -dependent error diffusion parameters */ 

TOKEMLUT *tokenLUT, /* Token LUT used in dot size diffusion */ 

short *dotdensityLUT, /* Optimal dot density curve */ 

short *dotsizeLUT) /* Optimal dot size curve */ 



{ 



short *fm_err , *am_err ; 

unsigned int i, j , token_img_width, mod_height; 

/* Initialize first row of fm error buffer */ 
srand(l) ; /* fix the seed */ 

fm_err = (short*) malloc (sizeof (short ) * (width+2 ) ) ; 
for(j = 0; j<width+2; j++) 

fm_err[j] = (rand () %128-64) ; /* Initialization */ 



/* initialize first row of am (dot size) error buffer */ 



am_err = (short*) malloc (sizeof (short) * ( (width>>l ) +4 ) ) ; 
for(i = 0; i < (width>>l ) +4 ; i++) 

am_err[i]=0; /* Initialization */ 

/* Avoid boundary because of pairwise row process */ 
if (height & 1) 

mod_height - height - 1; 
else 

mod_height = height; 

/* Process the input image with 2 rows each time */ 
for(i=0; i<mod_height ; i+=2) { 

if((i%600) == 0) printf ("amfm_ed: starting row %d\n" , i) ; 

amfm (width, i , contone_img [i] , token_img [i] , fm_err , am_err,dbs_screen, \ 
tdedpara, tokenLUT, dotdensityLUT, dotsizeLUT) ,- 

} 

/* Take care the last row if necessary */ 

if (height & 1) 

{ 

token_img_width = width/ 4; 
for ( j = 0 ; j <token_img_width; j ++) 

* (token__img [mod_height] + j ) = 0 ; 

} 



free (fm_err) ; 
free (am_err) ; 
return; 



} 

/* Define macro of FM_EVEN__ROW which does fm for a pair pixels in even row 
# define FM_EVEN_ROW \ 

/* First process FM (dot density) for left pixel in pixel pair. */ 

\ 

/* Get first pixel */ \ 
pixela = * (img_injptr++) ; \ 

\ 

/* Use look-up-table to get dot density */ \ 
dotdensity = dotdensityLUT [pixela] ; \ 

\ 

/* Compute look-up table entries for tone dependent error diffusion */ 
tded_ptr = (short*) (tdedpara + dotdensity) ; \ 
T2 = tded__ptr[0] 



DT = tded_ptr[l] 

Wl = tded_ptr[2] 

W2 = tded_ptr[3] 

W3 = tded_ptr[4] 

W4 = tded_ptr[5] 



/* compute dotdensity modified by diffused error */ 

mod_input = dotdensity + *f m_err jptr ; \ 

\ 

/* Threshold modifed dotdensity */ \ 
thresholding = mod_input - (dbs_jpat_rowptr [ j %SCREENWIDTH] * DT + T2); 
output = (thresholding > 0) ? 255 : 0; \ 

\ 



/* Compute weighted errors */ \ 
error = output - mod_input; 

el = (Wl * error) >>8; \ 

e2 = {W2 * error) >>8; \ 

e3 = (W3 * error) >>8; V 

/*e4 = (W4 * error) >>8 ; */ \ 

e4 = error - el - e2 -e3 ; \ 
/* Diffuse error forward in 1-D error buffer */ 

* ( - -f m_err__ptr ) -= e4 / \ 

* (++fm_err_jptr) = fm_tmp - e3 ; \ 
* (++fm_errj)tr) -= el; \ 
f m_tmp = -e2 ; \ 

/* Now process FM (dot density) for right pixel in pixel pair. */ 
/* Use same TDED parameters as for left pixel. */ 

\ 

/* Get second pixel */ \ 
pixelb = * (img_in_ptr++) ; \ 

/* Use look-up-table to get dot density */ 
dotdensity = dotdensityLUT [pixelb] ; 

\ 

mod_input = dotdensity + *fm_err_ptr; \ 
error = - mod input; /* suppress dot firing at this pixel */ 

\ 

el = (Wl * error) >>8; \ 

e2 = (W2 * error) >>8; \ 

e3 = (W3 * error) >>8; \ 

/*e4 = (W4 * error) >>8;*/ \ 

e4 = error - el - e2 -e3; \ 
/* Using the tded weights of the left pixel */ 

* ( - - fm_err_ptr) -= e4 ; \ 

* (++fm_err_j?tr) = f m_tmp - e3 ; \ 
* (++fra_err_ptr) -= el; \ 
f m_tmp - -e2; \ 

j += 2; 

/* Define macro of AM_ERROR_EVEN_ROW which computes and distributes 
dot size error for a pair pixels in even row */ 

#define AM_ERROR_EVEN_ROW \ 

/* el = Fl * error; */ \ 

e2 = F2 * error; \ 

e3 = F3 * error; \ 

e4 = (F4 * error) ; \ 

el = error*16 - e2 -e3 -e4; \ 

am_err__ptr - = 2; \ 

* (am__err j)tr) - =e4 ; \ 
* (++am_err_ptr) -=e3; \ 

* (++am__err_jptr) = -e2; \ 

* (++am__err_ptr) -=el; 

/* Define macro of FM_ODD_ROW which does fm for a pair pixels in odd 
#define FM_ODD_ROW \ 
/* First process FM (dot density) for right pixel in pixel pair */ 

\ 

/* Get right pixel */ \ 
pixela = * (img_in_ptr- - ) ; \ 



\ 

/* Use look-up-table to get dot density */ \ 
dotdensity = dotdensityLUT [pixela] ; \ 

\ 

/* Compute look-up table entries for tone dependent error diffusion */ \ 
tded_ptr = (short*) (tdedpara + dotdensity) ; \ 

T2 = tded_ptr [0] ; \ 

DT = tded_j>tr [1] ; \ 

Wl = tded__ptr [2] ; \ 

W2 = tded_j>tr [3] ; \ 

W3 = tded_jptr [4] ; \ 

W4 = tded_j>tr [5] ; \ 

/* Compute dotdensity modified by diffused error */ \ 
mod input = dotdensity + *fm_err_ptr; \ 

\ 

/* suppress this dot and compute the error */ \ 
error = - mod input ; \ 

\ 

/* Compute weighted errors */ \ 
el = (Wl * error) >>8; \ 
e2 = (W2 * error) >>8; \ 
e3 = (W3 * error) >>8; \ 
/*e4 = ( (W4 * error) >>8) ; */ \ 
e4 = error - el - e2 -e3; \ 

\ 

/* Diffuse error forward in 1-D error buffer */ \ 
* (++fm_err_ptr) - = e4; \ 
* (--fm_err_ptr) = fm_tmp - e3; \ 
* (--fm_err__ptr) el; \ 
fm tmp = -e2; \ 

\ 

/* Now process FM (dot density) for Left pixel in a pair */ \ 

\ 

/* Get second pixel */ \ 
pixelb = * (img_in_ptr- - ) ; \ 

/* Use look-up-table to get dot density */ \ 
dotdensity = dotdensityLUT [pixelb] ; \ 

\ 

mod input = dotdensity + *fm_err_ptr; \ 

\ 

/* Threshold modifed dot density */ \ 
thresholding = mod_input - (dbs_pat_rowptr [ ( j -1) %SCREENWIDTH] * DT + T2) ; \ 
output « (thresholding >0)?255:0; \ 

\ 

j -= 2; \ 

error = output - mod_input ; \ 

\ 

el = (Wl * error) >>8; \ 
e2 = (W2 * error) >>8; \ 
e3 = (W3 * error) >>8; \ 
/*e4 = (W4 * error) >>8;*/ \ 
e4 - error - el - e2 -e3; \ 

\ 



* (++fm_err_ptr) -= e4 ; \ 

* ( - -fm_err_ptr) = fm_tmp - e3 ; \ 

* ( - -fm_err _jptr) - = el; \ 
f m_tmp = -e2; 

/* Define macro of AM_ERROR_0DD_ROW which computes and distributes 

dot size error for a pair pixels in odd row */ 

#define AM_ERROR_ODD_ROW \ 

/* el = Fl * error; */ \ 

e2 = F2 * error; \ 

e3 = F3 * error; \ 

e4 = F4 * error; \ 

el = error*16 - e2 -e3 -e4; \ 

\ 

am__errjptr +=2; \ 

*am_errj)tr -= e4; \ 

* (--am_err_ptr) -= e3; \ 

* ( - -am_err_ptr) = -e2 ; \ 

* ( - -am_err_j)tr) - = el; 



/* This subroutine only processes 2 rows */ 
/* Assume width of image is multiple of 8 */ 
void amfm ( 

unsigned int width, /* Input image width */ 
unsigned int i, /* ith row */ 

unsigned char *img_in, /* ith row of input image array */ 

unsigned char *img_out, /* ith row of output image array */ 

short *fm_err, /* FM error buffer */ 

short *am_err, /* AM error buffer */ 

unsigned char ** dbs_screen, /* dbs_screen [SCREENHEIGHT] [SCREEWWIDTH] 

TD ED PARA *tdedpara, /* Tone -dependent error diffusion parameters */ 
TOKEMLUT *tokenLUT, /* Token look-up- table used in dot size diffusion */ 
short *dotdensityLUT, /* Optimal dot density curve */ 
short *dotsizeLUT) /* Optimal dot size curve */ 

{ 

short f m_tmp , thresholding ; 

short *fm_err_j?tr, *am_errj)tr; 

short pixela, pixelb, mod_dotsize, output; 

unsigned int j , img_out_width; 

unsigned char bit_pack; 

unsigned char *img_in_ptr , *img_out_ptr , *dbs_pat_rowptr ,- 
short dotdensity, dotsize, mod_input, error; 
short Wl, W2, W3 f W4 , T2 , DT, el, e2 , e3, e4 ; 
short *tded_ptr, *token_lut_ptr ; 
FILE *fp; 



/* */ 

/* serpentine even rows */ 
/* */ 

/* Initial points */ 
f m__tmp = 0 ; 

f m_err_jptr = fm_err+l ; 
am_errj)tr = am_err+2; 
img_in_ptr - img_in; 
img__out__ptr = img_out; 



/* Get row pointer of dbs pattern */ 

/* SCREENHEIGHT = 128, module 1 ( i + + ) %12 8 * can be replace by l (i + + )&127' */ 
dbsjpat_rowptr = dbs_screen [ (i++) %SCREENHEIGHT] ; 

/* Index through pixels in pairs */ 
for(j = 0; j<width;) { 

/* FM halftoning for first pixel pair */ 
FM_EVEN_ROW 

/* Begin section on AM halftoning for first pixel pair. */ 

/* This section computes 2-bit PWM codes for dot pairs. */ 

/* Errors in AM component are diffused using Floyd-Steinberg weights. */ 

/* Tokens of left and right pixels along with size error are precomputed */ 

/* and stored in tokenLUT */ 

/* Get diffused error from dot size error buffer */ 

/* This operation can be replace by y=x/16 without affecting */ 

/* the quality of the halftone. */ 

mod_dotsize = (*am_errj)tr+8) »4; 

bit_pack = 0; 
if (output) { 

mod__dotsize += dotsizeLUT [pixela] ; 

/* Get 2 PWM tokens and error corresponding to mod__dotsize */ 
/* Then pack the tokens */ 

token_lut_ptr = (short *) (tokenLUT + mod_dotsize) ; 

bit_pack = token__lut_ptr [0] <<6; /* Get and pack token for left pixel 

bit__pack += token_lut_ptr [1] <<4; /* Get and pack token for right pixel 
error = token_lut_ptr [2] ; /* Get size error for the pixel pair 

/ 

} 

else 

error = - mod__dotsize ; 

/* Compute and distribute dot size error */ 
AM_ERROR__EVEN_ROW 

/* FM halftoning for second pixel pair */ 
FM__EVEN_ROW 

/* Begin section on AM halftoning for second pixel pair. */ 
/* Same comments for AM halftoning of first pixel pair */ 
mod_dotsize = (*am_errj)tr+8) »4; 
if (output) { 

mod_dotsize += dotsizeLUT [pixela] ; 

/* Get 2 PWM tokens and error corresponding to mod_dotsize */ 
/* Then pack the tokens */ 

token_lut_ptr = (short *) (tokenLUT + mod_dotsize) ; 

bit_pack += token_lut_ptr [0] <<2; /* Get and pack token for left pixel 
bit_pack += token_lutjptr [1] ; /* Get and pack token for right pixel */ 



error = token_lut_ptr [2 ] ; 



/* Get size error for the pixel pair 



} 

else 

error = - mod_dotsize; 

/* Compute and distribute dot size error */ 
AM_ERROR_EVEN_ROW 

/* write the packed tokens to output image array */ 
* (img_out_ptr++) = bit_jpack; 

} /* end of ith row */ 

*/ 

serpentine odd rows */ 
*/ 

img_out_width = width/4; 
f m_tmp = 0 ; 

/* Set fm error buffer pointer to the end of fm_err buffer */ 

fm__err__ptr = fm_err+ width - 1; /* Offset by 1 */ 

/* Set am error buffer pointer to the end of am_err buffer */ 

am_err_ptr = am_err + (width>>l) ; /* Offset by 2 */ 

img_in_ptr = img_in+width*2 -2 ; 

img_out_ptr = img__out + img_out__width*2 - 1 ; 

/* Get row pointer of dbs pattern */ 

/* SCREENHEIGHT = 128, module ! i%128' can be replace by ' i&127 1 */ 
dbs j>at_rowptr = dbs_screen [i%SCREENHEIGHT] ; 

/* Index through pixels in pairs */ 

bit _pack = 0 ; 

for(j = width-2; j>2;) { 

/* FM halftoning for first pixel pair */ 
FM_ODD_ROW 

/* Begin section on AM halftoning for the first pixel pairl */ 

/* This section computes 2 -bit PWM codes for dot pairs. */ 

/* Errors in AM component are diffused using Floyd-Steinberg weights. */ 

/* Tokens of left and right pixels along with size error are precomputed */ 

/* and stored in tokenLUT */ 

/* Get diffused error from dot size error buffer */ 

/* This operation can be replace by y=x/16 without affecting */ 

/* the quality of the halftone. */ 

mod__dotsize = (*am_errjptr+8) »4; 

if (output) 

{ 

mod_dotsize += dot sizeLUT [pixelb] ; 

/* G e t 2 PWM tokens and error corresponding to mod_dotsize */ 
/* Then pack the tokens */ 

token_lut_ptr = (short *) (tokenLUT + mod_dotsize) ; 

error = token_lut_ptr [2 ] ; /* Get size error for the pixel pair */ 

bit_pack += token_lut_ptr [1] <<2 ; /* Get and pack token for right pixel 



bitjpack += token_lut_ptr [0] <<4 ; /* Get and pack token for left pixe 

} 

else 

error = - mod_dotsize; 
/* Compute and distribute dot size error */ 
AM_ERROR_ODDJROW 

/* FM halftoning for second pixel pair */ 
FM_ODD_ROW 

/* Section on AM halftoning for second pixel pair */ 

/* Same comments as on AM halftong for the first pixel pair */ 

mod_dotsize = (*am_errj)tr+8) >>4; 

if (output) 

{ 

mod_dotsize += dotsizeLUT [pixelb] ; 

/* Get 2 PWM tokens and error corresponding to mod_dotsize */ 
/* Then pack it */ 

token_lut_ptr = (short *) (tokenLUT + mod_dotsize) ; 

error = token_lut_ptr [2 ] ; /* Get size error for the pixel pair */ 

bit_jpack += token_lut__ptr [1] <<6 ; /* Get and pack token for right pixel 

/* Write the packed tokens to the output image array */ 
* ( img_out_pt r - - ) = bit_pack; 

bit_pack = token_lut__ptr [0] ; /* Get and pack token for left pixel */ 

} 

else 

{ 

* ( img_out_ptr- - ) = bit_pack; 

bit^pack = 0; 

error = - mod_dotsize; 

} 

/* Compute and distribute dot size error */ 
AM_ERROR_ODD_ROW 

} 

/* Take care of the most left three pixels of odd rows */ 

/* FM halftoning for the first pixel pair */ 
FM_ODD_ROW 

/* AM halftoning for first pixel pair */ 
mod_dotsize = (*am_errj>tr+8) >>4; 
if (output) 

{ 

mod_dotsize += dotsizeLUT [pixelb] ; 

/* Get 2 PWM tokens and error corresponding to mod_dotsize */ 
/* Then pack the tokens */ 

token__lutj>tr = (short *) (tokenLUT + mod_dotsize) / 
error = token_lut_ptr [2 ] ; /* Get size error for the pixel pair */ 
bit_pack += token_lut_ptr [1] <<2 ; /* Get token for right pixel */ 

bit_j>ack += token_lut_ptr [0] <<4 ; /* Get token for left pixel */ 

} 

else 

error - - mod_dotsize; 



/* Write the packed tokens to the output image array */ 
* ( img_out_ptr- - ) = bit_pack; 

/* Compute and distribute dot size error */ 
AM_ERROR_ODD_ROW 

return; 

} 

void get_stripe ( 
int img_height, 
int img_width, 

unsigned char **contone_img, 
int start j>oint, 
int stripe_width, 
unsigned char **stripe) 

{ 

int 

f or (i=0 ; i<img_height ; i++) 
for ( j = 0 ; j <stripe_width; j ++) 

stripe [i] [j] = contone_img [i] [ j +start_jpoint ] ; 

} 

void cut_out_result ( 
int img_height, 
int cut_offset, 
int cut__width, 
int store_of f set , 
unsigned char **output_stripe , 
unsigned char **output_img) 

{ 

int i , j ; 

for (i=0 ; i<img_height ; i++) 
f or ( j=0 ; j <cut__width; j + + ) 

output_img [i] [ j +store__of f set] =output_stripe [i] [ j +cut_of f set] 

} 



Appendix B2 



COMBINED DOT DENSITY AND DOT SIZE MODULATION 
Zhen He 

Charles A Bouraan 
Qian Lin 

10003284 
M-8658 US 



/* amf m_convert . c file */ 

/* 2 bits/pixel am/fm halftoning algorithm: part B */ 
/* Input a tiff file containing 2 -bit tokens */ 

/* Output a tiff file containing pulse width modulation codes */ 
/* Every pixel is left justified */ 

#include <stdio.h> 
^include <stdlib.h> 
#include <math.h> 
#include "tif f .h" 
#include "allocate . h" 
#define NE WRIGHT OxcO 
#define NEWLEFT 0x40 
#define NEWCENTER 0x00 

void token2pwm{int, int, unsigned char **, unsigned char **, unsigned char *) 
int main {int argc, char ** argv) 

{ 

unsigned char pwm [4] = { 0 , 21 , 42 , 63 } ; 

int i , j ; 
FILE * fp; 

struct TIFF_img input_img, output_img; 

if (argc<3 ) { 

printf ( "usage : %s token.tif output_img . tif \n" , argv [0] ) ; 
return 1 ; 

} 

printf ( "Mapping tokens to pulse width modulation codes. \n" ); 

/* read the input image */ 

if ( (fp=f open {argv [1] , "rb" ) ) ==NULL) { 

printf ("can not open file %s l\n" , argv [1] ) ; 

exit {2) ; 

} 

if (read_TIFF (fp r &input_img) ) { 

printf { "error reading input file\n n ); 
exit (3) ; 



} 

f close ( f p) ; 



/* Allocate memory for entire fm output image. */ 
get_TIFF{ &output_img , input_img - height , input_img . width* 4 , 'g' ); 

token2pwm ( input_img . he ight , input_img * width , input_img . mono , output_img . mono , 
pwm) ; 

/* write PWM codes image */ 

if{ (fp = fopen(argv[2] , "wb" ) ) ==NULL) { 

printf ("cannot open file %s\n" , argv[3]); 

exit (5) ; 

} 

if (write_TIFF (fp, &output_img) ) { 

printf ("\nError writing TIFF file %s\n", argv[2]); 
return 1; 

} 

f close (fp) ; 

/* free the space */ 

f ree_TIFF (& ( input_img) ) ; 

free_TIFF (& (output_img) ) ; 

f flush (stdout) ; 

return 0 ; 

} 



/* Map tokens to pulse width modulation codes */ 

void token2pwm ( 

int img__height, 

int img_width, 

unsigned char ** token, 

unsigned char ** output_pwm, 

unsigned char * pwm 

) 

{ 

int i , j , k ; 

for (i=0 ; i<img_height ; i++) 

for ( j = 0 ; j <img_width; j ++) 
{ 

output__pwm[i] [j*4] = pwm [ (token [i] [j ] &192) /64] + NEWLEFT; 
output_pwm[i] [j*4 + l] = pwm [ (token [i] [j]&48)/16] + NEWLEFT; 
output_pwm[i] [j*4+2] = pwm[ (token [i] [j]&12)/4] + NEWLEFT; 
output_pwm [i] [j*4 + 3] = pwm [token [i] [j]&3] + NEWLEFT ; 

} 

} 
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/* coef.h file */ 

#define S CREENHE I GHT 12 8 /* DBS screen height */ 

#define SCREENWIDTH 12 8 /* DBS screen width */ 

tdefine OVERLAP_W I DTH 16 /* Width of overlapping region 

#define Fl 0x0007 /* Floyd- Steinberg Weights 1/16 in Q4 */ 
#define F2 0x0003 /* Floyd- Steinberg Weights 5/16 in Q4 */ 
#define F3 0x0005 /* Floyd- Steinberg Weights 3/16 in Q4 */ 
#define F4 0x0001 /* Floyd-Steinberg Weights 7/16 in Q4 */ 
typedef struct TDEDPARA 

{ 

short T2; 
short DT; 
short Wl; 
short W2; 
short W3; 
short W4; 
} TDEDPARA; 

typedef struct TOKENLUT 
{ 

short left; 
short right; 
short s i z e_er ror ; 
} TOKENLUT; 

Static TDEDPARA TDEDcoef f [2 56 ] = { 



{76, 


0, 


181, 


0, 


3, 


72} , 


{76, 




181, 


0, 


3, 


72}, 


{79, 


o, 


172, 


1/ 


2 , 


81}, 


{80, 


0, 


161, 


14, 


18, 63}, 


{82, 


0, 


159, 


1, 


37, 




{83, 


0, 


149, 


6, 


5, 


96}, 


{83, 


o, 


141, 


30 , 


0, 


85}, 


{85, 


o, 


138, 


13, 


0, 


105} , 


{86, 


0, 


144 , 


10, 


1, 


101} , 


{85, 


0, 


129, 


48, 


3, 


76}, 


{86, 


0, 


123, 


31; 


1, 


101} , 


{87, 


o, 


123, 


29, 


3, 


101} , 


{87, 


0, 


115, 


28, 


5, 


108} , 



89, 0, 138, 19, 18, 81}, 
89, 0, 111, 17, 51, 77}, 
88, 0, 115, 31, 0, 110}, 

87, 0, 120, 16, 16, 104}, 

88, 0, 139, 12, 0, 105}, 

89, 0, 122, 19, 17, 98}, 

90, 0, 112, 32, 0, 112}, 

91, 0, 98, 34, 20, 104}, 
90, 10, 123, 16, 26, 9l}, 
93, 8, 126, 1, 74, 55}, 
89, 10, 89, 26, 71, 70}, 
89, 10, 89, 22, 43, 102}, 
89, 12, 91, 21, 34, 110}, 
88, 12, 85, 24, 30, 117}, 
88, 14, 85, 23, 30, 118}, 
84, 24, 113, 27, 13, 103}, 

82, 26, 113, 33, 0, 110}, 

83, 26, 109, 29, 9, 109}, 

84, 28, 106, 21, 29, 100}, 

85, 28, 103, 13, 56, 84}, 
96, 2, 102, 16, 57, 81}, 
93, 6, 102, 25, 28, 10l}, 

[91, 12, 102, 24, 32, 98}, 
>6, 2, 103, 24, 23, 106}, 
■94, 10, 99, 17, 62, 78}, 
[95, 6, 110, 12, 110, 24}, 
[97, 4, 114, 12, 112, 18}, 
[97, 6, 114, 11, 113, 18}, 
[96, 8, 111, 14, 110, 21}, 
[94, 12, 102, 17, 109, 28}, 
[94, 8, 79, 32, 108, 37}, 
[95, 6, 74, 35, 110, 37}, 
(97, 2, 70, 35, 111, 40}, 
(97, 4, 68, 33, 112, 43}, 
(97, 6, 69, 28, 112, 47}, 
(98, 6, 70, 22, 114, 50}, 
(97, 6, 68, 43, 113, 32}, 
(100, 4, 68, 22, 114, 52}, 
(99, 6, 71, 24, 112, 49}, 
{102, 2, 70, 23, 114, 49}, 
(100, 6, 68, 23, 114, 5l}, 
{100, 8, 66, 22, 116, 52}, 
{100, 8, 66, 24, 116, 50}, 
{96, 16, 75, 0, 122, 59}, 
{95, 16, 63, 0, 127, 66}, 
{95, 16, 56, 0, 130, 70}, 
{97, 14, 56, 0, 132, 68}, 
{97, 16, 59, 0, 132, 65}, 
{97, 16, 60, 0, 133, 63}, 
{98, 16, 62, 0, 133, 61} , 
{95, 26, 98, 0, 109, 49}, 
{97, 20, 65, 0, 132, 59}, 
{98, 18, 61, 0, 132, 63}, 
{99, 18, 63, 0, 131, 62}, 
{100, 16, 58, 0, 133, 65}, 
{100, 16, 58, 0, 131, 67}, 
{101, 16, 60, 0, 131, 65}, 



{101, 16, 63, 0, 129; 64}, 
{101, 16, 58, 0, 129, 69}, 
{102, 16, 71, 0, 123, 62}, 
{103, 8, 68, 23, 114, 51} , 
{103, 8, 66, 22, 116, 52}, 
{105, 6, 68, 22, 115, 5l}, 
{106, 4, 70, 22, 114, 50}, 
{108, 2, 69, 23, 113, 51} , 
{105, 8, 68, 22, 114, 52}, 
{108, 6, 70, 20, 115, 51}, 
{106, 8, 69, 27, 112, 48}, 
{109, 2, 65, 35, 112, 44}, 
{110, 4, 69, 34, 111, 42}, 
{110, 6, 72, 35, 110, 39}, 
{114, 0, 73, 34, 111, 38}, 
{110, 12, 94, 21, 108, 33}, 
{111, 12, 102, 15, 110, 29}, 
{116, 6, 114, 10, 113, 19}, 
{96, 16, 92, 16, 67, 81} , 
{100, 12, 95, 17, 67, 77}, 
{101, 12, 97, 19, 67, 73}, 
{99, 4, 101, 20, 45, 90}, 
{93, 4, 103, 25, 25, 103}, 
{94, 8, 101, 25, 33, 97}, 
{78, 24, 99, 26, 19, 112}, 
{81, 26, 104, 22, 24, 106}, 
{82, 26, 102, 26, 25, 103}, 
{91, 26, 109, 14, 46, 87}, 
{104, 10, 82, 0, 95, 79}, 
{107, 8, 83, 0, 97, 76}, 
{105, 8, 87, 2, 84, 83}, 
{81, 14, 86, 27, 25, 118}, 
{99, 12, 122, 0, 37, 97}, 
{102, 10, 117, 0, 45, 94}, 
{103, 10, 90, 21, 64, 81}, 
{105, 12, 122, 4, 51, 79}, 
{101, 12, 126, 9, 29, 92}, 
{88, 12, 121, 25, 0, 110}, 
{85, 12, 114, 25, 1, 116}, 
{89, 10, 109, 23, 10, 114}, 
{86, 12, 112, 29, 1, 114}, 
{89, 12, 119, 31, 0, 106}, 
{94, 10, 123, 37, 1, 95}, 
{93, 8, 117, 63, 1, 75}, 
{99, 6, 118, 75, 9, 54}, 
{97, 6, 120, 43, 3, 90}, 
{111, 6, 121, 35, 32, 68}, 
{95, 6, 116, 54, 0, 86}, 
{107, 6, 125, 39, 15, 77}, 
{93, 34, 137, 27, 19, 73}, 
{85, 44, 139, 33, 16, 68}, 
{87, 48, 146, 31, 23, 56}, 
{87, 44, 148, 22, 10, 76}, 
{93, 40, 152, 22, 11, 71} , 
{97, 44, 159, 4, 28, 65}, 
{95, 42, 161, 25, 4, 66}, 
{103, 48, 176, 3, 44, 33}, 



{101, 56, 165, 27, 55, 9}, 
{97, 56, 165, 27, 55, 9}, 
{103, 48, 176, 3, 44, 33}, 
{117, 42, 161, 25, 4, 66}, 
{113, 44, 159, 4, 28, 65}, 
{121, 40, 152, 22, 11, 71}, 
{123, 44, 148, 22, 10, 76}, 
{119, 48, 146, 31, 23, 56}, 
{125, 44, 139, 33, 16, 68}, 
{127, 34, 137, 27, 19, 73}, 
{141, 6, 125, 39, 15, 77}, 
{153, 6, 116, 54, 0, 86}, 
{137, 6, 121, 35, 32, 68}, 
{151, 6, 120, 43, 3, 90}, 
{149, 6, 118, 75, 9, 54}, 
{153, 8, 117, 63, 1, 75}, 
{150, 10, 123, 37, 1, 95}, 
{153, 12, 119, 31, 0, 106}, 
{156, 12, 112, 29, 1, 114}, 
{155, 10, 109, 23, 10, 114}, 
{157, 12, 114, 25, 1, 116}, 
{154, 12, 121, 25, 0, 110}, 
{141, 12, 126, 9, 29, 92}, 
{137, 12, 122, 4, 51, 79}, 
{141, 10, 90, 21, 64, 81}, 
{142, 10, 117, 0, 45, 94}, 
{143, 12, 122, 0, 37, 97}, 
{159, 14, 86, 27, 25, 118}, 
{141, 8, 87, 2, 84, 83}, 
{139, 8, 83, 0, 97, 76}, 
{140, 10, 82, 0, 95, 79}, 
{137, 26, 109, 14, 46, 87}, 
{146, 26, 102, 26, 25, 103}, 
{147, 26, 104, 22, 24, 106}, 
{152, 24, 99, 26, 19, 112}, 
{152, 8, 101, 25, 33, 97}, 
{157, 4, 103, 25, 25, 103}, 
{151, 4, 101, 20, 45, 90}, 
{141, 12, 97, 19, 67, 73}, 
{142, 12, 95, 17, 67, 77}, 
{142, 16, 92, 16, 67, 8l}, 
{132, 6, 114, 10, 113, 19}, 
{131, 12, 102, 15, 110, 29}, 



{132 , 


12 , 


94 , 


21, 


, 108, 


33} 


{140, 


0, 


73 , 


34 , 


111, 


38}, 


{138, 


6, 


72, 


35, 


110, 


39}, 


{140, 


4 , 


69, 


34, 


111, 


42}, 


{143 , 


2, 


65, 


35, 


112 , 


44}, 


{140 , 


8, 


69, 


27 , 


112 , 


48} , 


{140 , 


6 , 


70 , 


20, 


115, 


51}, 


{141, 


8, 


68 , 


22 , 


114, 


52}, 


{144 , 


2 , 


69, 


23, 


113, 


51}, 


{144 , 


4 , 


70 , 


22, 


114 , 


50}, 


{143, 


6 , 


68, 


22 , 


115, 


51}, 


{143, 


8, 


66 , 


22 , 


116, 


52} , 


{143, 


8, 


68 , 


23 , 


114, 


51}, 


{136, 


16, 


71 , 


0 , 


123, 


62} , 



{137, 16, 58, O f 129, 69}, 
{137, 16, 63, 0, 129, 64}, 
{137, 16, 60, 0, 131, 65}, 
{138, 16, 58, 0, 131, 67}, 
{138, 16, 58, 0, 133, 65}, 
{137, 18, 63, 0, 131, 62}, 
{138, 18, 61, 0, 132, 63}, 
{137, 20, 65, 0, 132, 59}, 
{133, 26, 98, 0, 109, 49}, 
{140, 16, 62, 0, 133, 6l}, 
{141, 16, 60, 0, 133, 63}, 
{141, 16, 59, 0, 132, 65}, 
{143, 14, 56, 0, 132, 68}, 
{143, 16, 56, 0, 130, 70}, 
{143, 16, 63, 0, 127, 66}, 
{142, 16, 75, 0, 122, 59}, 
{146, 8, 66, 24, 116, 50}, 
{146, 8, 66, 22, 116, 52}, 
{148, 6, 68, 23, 114, 5l}, 
{150, 2, 70, 23, 114, 49}, 
{149, 6, 71, 24, 112, 49}, 
{150, 4, 68, 22, 114, 52}, 
{151, 6, 68, 43, 113, 32}, 
{150, 6, 70, 22, 114, 50}, 
{151, 6, 69, 28, 112, 47}, 
{153, 4, 68, 33, 112, 43}, 
{155, 2, 70, 35, 111, 40}, 
{153, 6, 74, 35, 110, 37}, 
{152, 8, 79, 32, 108, 37}, 
{148, 12, 102, 17, 109, 28}, 
{150, 8, 111, 14, 110, 21}, 
{151, 6, 114, 11, 113, 18}, 
{153, 4, 114, 12, 112, 18}, 
{153, 6, 110, 12, 110, 24}, 
{150, 10, 99, 17, 62, 78}, 
{156, 2, 103, 24, 23, 106}, 
{151, 12, 102, 24, 32, 98}, 
{155, 6, 102, 25, 28, 10l}, 
{156, 2, 102, 16, 57, 8l}, 
{141, 28, 103, 13, 56, 84}, 
{142, 28, 106, 21, 29, 100}, 
{145, 26, 109, 29, 9, 109}, 
{146, 26, 113, 33, 0, 110}, 
{146, 24, 113, 27, 13, 103}, 
{152, 14, 85, 23, 30, 118}, 
{154, 12, 85, 24, 30, 117}, 
{153, 12, 91, 21, 34, 110}, 
{155, 10, 89, 22, 43, 102}, 
{155, 10, 89, 26, 71, 70}, 
{153, 8, 126, 1, 74, 55}, 
{154, 10, 123, 16, 26, 9l}, 
{163, 0, 98, 34, 20, 104}, 
{164, 0, 112, 32, 0, 112}, 
{165, 0, 122, 19, 17, 98}, 
{166, 0, 139, 12, 0, 105}, 
{167, 0, 120, 16, 16, 104}, 
{166, 0, 115, 31, 0, 110}, 



{165, 0, 111, 17, 51, 77}, 

{165, 0, 138, 19, 18, 81}, 

{167, 0, 115, 28, 5, 108}, 

{167, 0, 123, 29, 3, 10l}, 

{168, 0, 123, 31, 1, 101}, 

{169, 0, 129, 48, 3, 76}, 

{168, 0, 144, 10, 1, 101}, 

{169, 0, 138, 13, 0, 105}, 

{171, 0, 141, 30, 0, 85}, 

{171, 0, 149, 6, 5, 96}, 

{172, 0, 159, 1, 37, 59}, 

{174, 0, 161, 14, 18, 63}, 

{175, 0, 172, 1, 2, 81}, 

{178, 0, 181, 0, 3, 72}, 

{178, 0, 181, 0, 3, 72}, 

}; 

static short OptSizeLUT [256] 

120, 

118, 

117, 

116, 

115, 

114, 

112, 

111, 

109, 

108, 

107, 

105, 

104, 

102, 

101, 

100, 

100, 

98, 

97, 

97, 

96, 

95, 

94, 

93, 

93, 

92, 

91, 

90, 

89, 

89, 

88, 

87, 

86, 

86, 

85, 

84, 

84, 

83, 

82, 

82, 



81, 
80, 
80, 
79, 
79, 
78, 
78, 
78, 
77, 
77, 
77, 
76, 
76, 
76, 
75, 
75, 
75, 
75, 
74, 



69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 



o 
m 
ru 

13 

o 



o 
ffi 

If! 
SI 
CI 

o 



74, 
74, 
73, 
73, 
72, 
72, 
72, 
71, 
71, 
71, 
70, 
70, 
70, 
70, 
70, 
70, 
70, 
69, 



69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 

O 69 ' 



69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 



o 

03 
Pi 
4= 
O 
Q 



ffl 
,p 

Ul 
M 
m 

a 



69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 



o 
m 
ui 

M 

€i 

o 

s 

D 
03 

-p 

O 

a 



69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

69, 

68, 

68, 

68, 

68, 

67, 

67, 

67, 

67, 

66, 

66, 

65, 

65, 

65, 

64, 

64, 

63, 

63, 

62, 

62, 

61, 

61, 

60, 

60, 

59, 

59, 

58, 

58, 

57, 

57, 

56, 

56, 

55, 

55, 



54, 
54, 
53, 
53, 
52, 
52, 
52, 
51, 
51, 
50, 
50, 
49, 
49, 
48, 
48, 
47, 
47, 
47, 
46, 

O 46 ' 

to 46 ' 
45, 



4= 

Ul 



44, 
44, 
44, 
43, 

O 43, 
« 43, 
p 42, 

42, 
42, 
41, 
41, 



ni 
a 

O 41, 
40, 



40, 
40, 
39, 
39, 
38, 
38, 
38, 
38, 

}; 

static short OptDensityLUT [256 ] 

128, 

125, 

123, 

121, 

119, 

117, 

116, 

115, 

114, 

113, 



113, 
112, 
112, 
111/ 
HI, 
HI, 
111, 
111, 
110, 
110, 
110, 
110, 
110, 
110, 
109, 
109, 
109, 
109, 
109, 

O 109 ' 

*3 109, 

m 108 - 

108, 



108, 
108, 
108, 



o 

O 



m 

108, 

O 108, 
~ 108, 

Q 107 ' 
m io7, 

ft 107 < 

107, 
107, 
107, 
107, 
107, 
107, 
107, 
107, 
107, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
106, 
105, 
105, 
105, 



Q 
m 

m 

UI 
Si 
« 

= 

O 

m 
m 

-c 

o 
a 



105, 

105, 

105, 

104, 

104 , 

104, 

104, 

104, 

103 , 

103, 

103, 

102, 

102 , 

102, 

101, 

101, 

101, 

100, 

100, 

100, 

99, 

99, 

98, 

98, 

98, 

97, 

97, 

96, 

96, 

96, 

95, 

95, 

94, 

94, 

93, 

93, 

93, 

92, 

92, 

91, 

91, 

91, 

90, 

90, 

89, 

89, 

88, 

88, 

88, 

87, 

87, 

86, 

86, 

85, 

85, 

84 , 

84, 



o 
m 

II! 
Si 

m 
a 

~ 

o 
m 

fU 

o 

a 



83, 

83, 

83, 

82 , 

82, 

81, 

81, 

80, 

80, 

79, 

79, 

78, 

78, 

77 , 

77, 

76, 

76, 

75, 

75, 

74, 

74, 

73, 

73, 

72 , 

71, 

71, 

70, 

70, 

69, 

69, 

68, 

68, 

67, 

67 , 

66, 

66, 

65, 

65, 

64, 

64, 

64, 

63, 

63, 

62 , 

62, 

62 , 

61, 

61, 

61, 

60, 

60, 

60, 

60 , 

59, 

59, 

59, 

59, 



o 
m 

4= 

m 
%J 
m 
o 

s 

a 
en 
ru 
4= 

a 
a 



59, 

58, 

58, 

58, 

58, 

58, 

58, 

58, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

57, 

56, 

56, 

56, 

56, 

56, 

56, 

56, 

56, 

56, 

56, 

55, 

55, 

55, 

55, 

55, 

55, 

54, 

54, 

54, 

54, 

53, 

53, 

53, 

52, 

52, 

52, 

51, 

51, 

50, 

50, 

49, 

48, 

48, 

47, 

46, 

45, 

44, 



43, 

42, 

41, 

40, 

39, 

37 , 

36, 

35, 

33, 

31, 

29, 

27, 

24, 

20, 

16, 

11, 

7, 

0, 

}; 

Static TOKENLUT TokenLUT [18 0 ] - { 





o. 


30}, 


/* 


-30 


*/ 


{o, 


o, 


29} , 


/* 


-29 


*/ 


{0, 


0, 


28} , 


/* 


-28 


*/ 


{o. 


0, 


27}, 


/* 


-27 




{o. 


0, 


26}, 


/* 


-26 


*/ 


{o, 


0, 


25} , 


/* 


-25 


*/ 


{o, 


0, 


24}, 


/* 


-24 


*/ 


{o. 


o. 


23}, 


/* 


-23 


*/ 


{o, 


o. 


22} , 


/* 


-22 


*/ 


{o. 


o. 


21}, 


/* 


-21 




{o. 


o. 


20} , 


/* 


-20 


*/ 


{o, 


0, 


19}, 


/* 


-19 




{o, 


o, 


18}, 


/* 


-18 


*/ 


{o. 


0, 


17}, 


/* 


-17 


*/ 


{o. 


o. 


16}, 


/* 


-16 


*/ 


{o. 


0, 


15}, 


/* 


-15 


*/ 


{o, 


o, 


14}, 


/* 


-14 


*/ 


{o. 


0, 


13}, 


/* 


-13 


*/ 


{o. 


o, 


12}, 


/* 


-12 




{o. 


0, 


11}, 


/* 


-11 




{o. 


0, 


10}, 


/* 


-10 




{o. 


0, 


9}, 


/* 


-9 




{o, 


0, 


8}, 


/* 


-8 




{o. 


o, 


7}, 


/* 


-7 


*/ 


{o. 


0, 


6}, 


/* 


-6 


*/ 


{o. 


0, 


5}, 


/* 


-5 


*/ 


{o, 


o, 


4}, 


/* 


-4 




{o, 


o, 


3}, 


/* 


-3 


*/ 


{o. 


o, 


2}, 


/* 


-2 


*/ 


{o, 


o, 


1}, 


/* 


-1 




{o, 


o, 


o}, 


/* 


0 * 


/ 



{0, 0, -1}, /* 1 */ 

{0, 0, -2}, /* 2 */ 

{0, 0, -3}, /* 3 */ 

{0, 0, -4}, /* 4 */ 

{0, 0, -5}, /* 5 */ 

{0, 0, -6}, /* 6 */ 



lo. 


o , 


-7), /* 


7 


*/ 






0 , 


-8}, /* 


8 


*/ 






o , 


-9}, /* 


9 


*/ 




(o, 


o , 


-10 } , 




/* 


10 


(0, 


0 , 


-11} , 




/* 


11 


f 0 


o , 


-12 } , 




/* 


12 


{o, 


o, 


-13}, 




/* 


13 


{o, 


o, 


-14}, 




/* 


14 


{o, 


o, 


-15}, 




/* 


15 


{o. 


o, 


-16} , 




/* 


16 


{o, 




-iv}, 




/* 


17 


{o, 




-18}, 




/* 


18 


{o. 




-19}, 




/* 


19 


{o, 


o, 


-20} , 




/* 


20 


{o. 




-21} , 




/* 


21 



\ z * 


o 


9 0 \ 




22 




I 9 


r\ 
U / 






23 


/ 


[ 9 
1 2 ' 


n 


xo j f 




24 


/ 


/ 9 


n 


1 7 \ 


/* 


2 5 


*/ 


! i 
\ 2 ' 


ri 

u , 


1 fi\ 

J- 1> J , 


/ 


2 6 


/ 


I 9 




XD j , 


/* 


2 7 


* / 


t 9 




14] 


/* 


2 8 




io 
\ Z ' 


n 


XJ / , 




2 9 


/ 


\ Z ' 


w / 


19^ 


i 


3 0 


/ 


f 9 


A 
v / 


XX ] , 




3 1 




l*> 
\ Z / 


r\ 

U / 


xu / , 




32 


*/ 


/ 9 


n 


y / ' 




33 


/ 


I 9 


n 


a* 




34 




to 


U / 


7) 




35 




\ z / 


0 


6) 




36 


*/ 






5\ 
~> / , 




37 


*/ 


f 2 


o 


** j 1 


/* 


38 


*/ 


(2 


o 


0 J r 


1* 


39 


*/ 


{2, 


o, 


2}, 


/* 


40 


*/ 


{2, 


o, 


1}, 


/* 


41 


*/ 


{2, 


o, 


0}, 


/* 


42 


*/ 


{2, 


o, 


-1}. 


/* 


43 


*/ 


{2, 


o, 


-2}, 


/* 


44 


*/ 


{2. 


o, 


-3}, 


/* 


45 


*/ 


{2, 


o, 


-4}, 


/* 


46 


*/ 


{2, 


o, 


-5}, 


/* 


47 


*/ 


{2, 


o, 


-6}, 


/* 


48 


*/ 


{2, 


o, 


-7}, 


/* 


49 


*/ 


{2, 


o, 


-8}, 


/* 


50 


*/ 


{2, 


o, 




/* 


51 


*/ 


{2, 


o, 


-10} 






/* 


{3- 


0, 


10}, 


'/* 


53 


*/ 


{3, 


o, 


9}. 


/* 


54 


*/ 


{3, 


o, 


8}, 


/* 


55 


*/ 


(3, 


o, 


7}, 


/* 


56 


*/ 


{3, 


0, 


6}, 


/* 


57 


*/ 


{3, 


o, 


5}, 


/* 


58 


*/ 


{3, 


o, 


4}, 


/* 


59 


*/ 


{3, 


o, 


3}, 


/* 


60 


*/ 


{3, 


o, 


2}, 


/* 


61 


*/ 


{3, 


0, 


1}. 


/* 


62 


*/ 


{3, 


o, 


0}, 


/* 


63 


*/ 



{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 0, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 1, 
{3, 2, 



-1}, /* 64 */ 
-2}, /* 65 */ 
-3}, /* 66 */ 
-4}, /* 67 */ 
-5}, /* 68 */ 
-6}, /* 69 */ 
-7}, /* 70 */ 
-8}, /* 71 */ 
-9}, /* 72 */ 
-10}, /* 73 */ 

74 */ 

75 */ 

76 */ 

77 */ 

78 */ 

79 */ 

80 */ 

81 */ 

82 */ 

83 */ 

84 */ 

85 */ 

86 */ 

87 */ 

88 */ 

89 */ 

90 */ 

91 */ 

92 */ 

93 */ 
/* 94 */ 

95 */ 



10} , 
9}, 
8}, 
7}, 
6}, 
5}, 
4}, 
3}, 
2}, 
1}, 
0}, 

-2}, 
-3}, 
-4}, 
-5}, 
-6}, 
-7}, 
-8}, 

-9}, 
-10}, 
10}, /* 



/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 



{3, 


2, 


9}, 


/* 


96 * 


/ 


{3, 


2, 


8}, 


/* 


97 * 


/ 


{3, 


2, 


7}, 


/* 


98 * 


7 


{3, 


2 , 


6}, 


/* 


99 * 


7 


{3, 


2, 


5}, 


/* 


100 


*/ 


{3, 


2, 


4}, 


/* 


101 


*/ 


{3, 


2, 


3}, 


/* 


102 


*/ 


{3, 


2, 


2}, 


/* 


103 


*/ 


{3, 


2, 


1}, 


/* 


104 


*/ 


{3, 


2, 


o}, 


/* 


105 


*/ 


{3, 


2, 


-1}- 


/* 


106 


*/ 


{3, 


2 , 


-2}, 


/* 


107 


*/ 


{3, 


2, 


-3}, 


/* 


108 


*/ 


{3, 


2, 


-4}, 


/* 


109 


*/ 


{3, 


2, 


-$}, 


/* 


110 


*/ 


{3, 


2, 


-6}, 


/* 


111 


*/ 


{3, 


2, 


-7}, 


/* 


112 


*/ 


{3, 


2, 


-8}, 


/* 


113 


*/ 


{3, 


2 , 


-9}, 


/* 


114 


*/ 


{3, 


2, 


-10} , 


/* 


{3, 


3, 


10} , 


/* 


116 


*/ 


{3, 


3 , 


9}- 


/* 


117 


*/ 


{3, 


3, 


8}, 


/* 


118 


*/ 


{3. 


3 , 


7}- 


/* 


119 


*/ 


{3, 


3 , 


6}, 


/* 


120 


*/ 





3 

° t 


5) 


/* 

f 


121 * / 






(3 


3 




f 


122 */ 






l 3 / 


3 


3 ) 




123 */ 






(3 
i ' 


3 






124 * / 






(2 
i ^ / 


3 






125 */ 






(3 


3 


0} 


/* 


126 */ 






(3 
\ J ' 


3 


J ' 


/* 


127 */ 






(3 


3 


-2) 




128 */ 






(3 


3 


-3) 


/* 


12 9 */ 






(3 


3 


-4}, 


/* 


130 */ 






(3 


3 


-5}, 


/* 


131 * / 






(3 


3 




/* 


132 */ 






(3 


3 / 


-7}, 


/* 


133 */ 






13 


3 

J r 


-8}, 


/* 


134 * / 






(3 


3 


-9}, 


/* 


135 */ 






(3 


3 t 


-10}, 






136 




(3 


3 


-11}. 




/* 


137 


*/ 


(3 


3 


-12}, 


/* 
/ 


13 8 


/ 


i 3 


3 


-13}, 




/ 


139 


* / 


1 3 


3 


-14}, 




/* 
/ 


140 




{3, 


3 , 


-15}, 




/* 


141 


*/ 


(3 


3 t 


-16}, 




/* 


142 


*/ 


{3, 


3, 


-17} 




/* 


143 


*/ 


{3, 


3, 


-18} 




/* 


144 


*/ 


(3, 


3, 


-19} 




/* 


14 5 


*/ 


{3, 


3; 


-20} 


/ 


/* 


146 


*/ 


{3, 


3 , 


-21} 




/* 


147 


*/ 


{3, 


3, 


-22 } 




/* 


148 


*/ 


{3, 

}; 


3 , 


-23} 




/* 


149 


*/ 



